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Abstract. Feynman using a formula, known as Feynman's law for a moving 
point charge, explained the phenomenon of synchrotron radiation and derived 
formulas for phenomena, concerning electromagnetic radiation at large dis- 
tance from the source, such as reflection, refraction, interference, diffraction, 
and scattering. These facts show the importance of the formula. 

The formula is supposed to represent the intensities of the electric and 
magnetic fields in free space, that is satisfying homogeneous system of Maxwell 
equations. 

Feynman's law contains a mathematical inconsistency. It involves implicitly 
a field representing the retarded time and an ordinary derivative in place where 
a partial derivative should be. 

In this note the author shows how to prove, using Banach's contraction 
mapping theorem, the existence and uniqueness of the retarded time field for 
any relativistically admissible trajectory of a point mass. 

The Lorentzian frame, the trajectory, and the retarded time field uniquely 
determine a system of fundamental fields. 

By means of these fields one can represent and establish relations between 
the system of wave-gauge equations and Maxwell equations, and to prove that 
electromagnetic field represented by the amended Feynman's formula satisfies 
the homogenous system of Maxwell equations and the system of wave-gauge 
equations. 

As applications the author proves the existence and uniqueness of the solu- 
tion to the n-body problem in the resulting joint gravitational and electromag- 
netic fields as also in any other relativistic force field representing a regular 
nonanticipating operator of the system of trajectories. 



1. Analysis of the original Feynman's formula 
for a moving point charge 

Importance of Maxwell's laws for logical analysis of electromagnetic phenomena 
stems from Einstein's work on special theory of relativity 8 . An asymmetry in the 
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interpretation of an eleetromagnetic phenomena based on Lorentz force [14j , acting 
on a charged particle, lead Einstein to the discovery of this theory. 

Einstein proved that laws of physics have to be invariant under transformations 
discovered earlier by Lorentz [Mj. In the process he established that Maxwell's 
equations are invariant under such transformations [SI. 

Feynman derived heuristically a formula for the electromagnetic field generated 
by a moving point charge. In this paper we will prove that after replacing the 
ordinary derivative with respect to time by the partial derivative and treating all 
quantities appearing in Feynman's formula as fields, the corrected formula yields 
electromagnetic field {E, B) satisfying Maxwell equations. In the process we shall 
also prove that Lienard-Wiechert formulas provide an electromagnetic potential for 
such a field. 

To see that one has to be careful when treating derivatives 

d , d 
— and — 
dt dt 

of composed function consider the following example. Let u = u{x^y,z,t) be a 
differentiable function in i?^ representing some physical quantity as a function of 
position and time. Let {x{t),y{t), z{t)) denote a path of a particle as a function of 
time. Then the derivative with respect to time along the path is 

d ... 

— U = Ur,. X + Uy y + U;, Z + Tij. 

Obviously it is not the same as partial derivative on the path 

d 



Any moving point mass mp obeying Einstein's laws of special theory of relativity 
generates a pair of fields {E, B) satisfying Maxwell's equations just from geometrical 
considerations of the Lorentzian frame. No other physical properties are required. 
Thus the obtained field E may as well represent for instance the illusive gravity 
field that Einstein was looking for [9|. 

To formulate briefiy the main result of the paper assume that in a given Lorentzian 
frame we introduce fundamental fields that are uniquely determined by the trajec- 
tory of a moving point mass. By means of these fields and their partial derivatives 
we can represent in explicit form fields given by amended Feynman's and Lienard- 
Wiechert's formulas and prove the relations between such fields and Maxwell's and 
wave equations. 

The main result can be stated as follows: 

Theorem 1.1 (Bogdan- Feynman Theorem). Assume that in a Lorentzian frame 
is given an admissible trajectory t^r2 (t) from R into representing a path of a 
point mass mp. Let G denote the set of points that do not lie on the path. Define 
Newton- Feynman field by the formula 

Id 19^ 
E{ri,t)^u^e+-u-^ — {u'^e) + — — {e) for all n € R^ , t € R 

and the associated field by 

B{ri,t)^-exE for all rieR^,teR. 
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Every pair of such fields {E, B) satisfies the homogenous system of Maxwell 
equations and the homogenous system of wave equations on the set G. As such 
these fields propagate through space with velocity c of light. 

The ^-vector obtained by Lienard- Wiechert formulas provides a pair of scalar- 
vector potentials for the pair {E, B). 



For the notion of an admissible trajectory see (|2.1|) and for definitions of the 
involved fields u and e see (14.11). 



We shall start with analysis of Feynman's original formulas: 

The intensity of the electric field E and of the magnetic field B at any time t 
and any point ri G R^, not lying on the trajectory of a moving charge q, are given 
by 



(1-1) E=-^ 
Alien 



e \r\ d 
rP c dt 





1 d2 - 




C2 df^ ^ 



1 

B ^ - exE, 



where r — ri — r2{t') is a vector, starting on the trajectory t i— > r2{t) of a moving 
point charge at the retarded time t' = t — |r|/c and ending at the point ri at time t 
where the field is to be evaluated. Here e denotes the unit vector corresponding to 
the vector r, and c the speed of light, and eg th^ electrostatic constant. 

For reference see Feynman-Leighton-Sands [10], vol. 1, chapter 28, formulas 
(28.3) and (28.4). 



Let us pose for a moment to analyze the formulas 11.11 for their mathematical 
content. On the left side of the formula we have the quantity E representing 
intensity of the electric field at point ri at time t. Thus E — E{ri,t) represents a 
function of the point (ri, t). So on the right side we should also find an expression 
representing a function of these variables. The vector r — ri — r2{t') and therefore 
the vector e — r/\r\ involves a retarded time that must satisfy the equation 

t' ^t^\r\/c 

equivalently 

t' = t-\r,-r2{t')\lc 

Since the trajectory t i— >■ r2{t) is fixed the above equation implicitly defines the 
retarded time as a function of (ri,t). If one can solve the above equation explicitly 
we will get t' = t'{ri,t) in the form of a function. Once this function is known 
and we can prove that it has continuous partial derivatives up to order 2 and the 
trajectory r2{t) itself has continuous derivative up to order 2, then the composed 
functions 

r{ri,t) = ri - r2{t' {ri,t)) and e{ri,t) ^ ——^——r{ri,t) 

\nri,t)\ 

will be well defined on the set G = {{ri,t) : \r{ri,t)\ > 0} and the functions will 
be of class C^. 



Our immediate goal will be to find an explicit formula for the retarded time 
function and to prove that the functions appearing in Feynman's formula indeed 
have the required properties. 
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To simplify the notation we select units of measure so that the speed of light is 
c — 1 and the electrostatic constant satisfies the condition Aireo — 1. We shall also 
assume that all the intensities of the involved fields are per unit of charge, that is 
we assume that g = 1. Let us introduce also function u = 

After these modifications and clarifications the formula for E can be reduced to 
a formula 

Definition 1.2 (Newton- Feynman formula). 

(1.2) E{ri,t) = u^e + u-^ — {u^e) + — e for all (ri,t)eG. 

Notice that, when the point mass representing the charge, is in its rest frame, 
that is r2(i) — for all t, then the above formula reduces to a single term 

1 



u e = -{n - ra) 

Fi - ^2 



differing from Newton's gravitation formula just by a constant of proportionality. 
Thus it is proper to call the expression (|1.2I) the Newton- Feynman formula. 



2. Considerations concerning trajectories 



In this section we shall refine the notion of the trajectory t i— > r2{t) so that 
we could prove that the functions involved in amended Feynman's formulas are 
of class that is the functions and their partial derivatives up to order 3 are 
continuous on the set G of points not lying on the trajectory. We will need this 
class of regularity in order to be able to prove that the fields generated by amended 
Feynman's formula coincide with fields generated by means of Lienard-Wiechert 
potentials. 

A trajectory of a path of a point mass can be parameterized in several different 
ways. It is important to understand which of these parameterizations depend on 
the Lorentzian frame, which are invariant under Lorentzian transformations and 
thus belong to Einstein's special theory of relativity, and which can be carried over 
to general theory of relativity. 

Following Einstein [B] we define a Lorentzian frame to consist of an orthogonal 
coordinate system in R^, having right hand orientation of axes, and equipped with 
a clock. A physical event in such a frame is described by a point (ri,i), where ri 
denotes a position in R'^ and t G R the time of the event. Denote such a frame by S. 
Assume that S' denotes another Lorentzian frame whose origin initially coincides 
with the origin of the frame S. Moreover the frame S' moves as a rigid body away 
from the frame of S* at a constant velocity. The transformation of coordinates 
of events from the frame S into the frame S' forms a linear transformation that 
preserves the quadratic form 

\r[\'-itr = \n\'~t'. 

By geometry of Lorentz space-time we shall understand the product space 
R^ X R with the transformations of coordinates as described above. These trans- 
formations form a group with composition of transformations as a group operation. 
Though one could expand the group by adding affine transformations, the linear 
transformations are sufficient for description of dynamics in physical processes in 
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Einstein's special theory of relativity. Any afBne orthogonal transformation can be 
reduced to a linear one just by moving the origin of the coordinate system. 

More general groups of transformations related to Lorentz group were studied by 
several authors. For generalization of such transformations and further references 
see Vogt [16]. 

Let a I— > y{a) be a mapping of an interval / into R^, of class C^, that is having 
continuous derivatives up to order 3 on the entire interval /. Assume that the 
mapping forms a parametric representation of a path of a point mass in space. 

Moreover assume that the tangent vector field y' consists of time-like vectors 
that is 

(2.1) y[{a)'^ + y'^iaf + y'^iaf - y'^ia)'^ < for all a e R. 

As a derivative of a covariant field with respect to a free parameter the tensor 
y'jia) itself forms a covariant field over /. Thus it is invariant under Lorentzian 
transformations and it can be carried over to the general theory of relativity as in 
Dirac |7] . The transition from covariant to contravariant tensors is given by means 
of transformation 

y^ = g'^'yk 

where summation is with respect to index k = 1,2,3,4 and the matrix g^^ for 
orthogonal axes has elements on the diagonal equal respectively 1,1,1,-1 and non- 
diagonal elements are zero. 

The time along the path is given by t = 2/4(0;), since from the relation (|2.ip 
follows that that ^ = y'^{a) > for all a £ I, the correspondence a 1-^ t represents 
a diffeomorphism of / onto some interval J and is also of class C^, that is both 
maps a t-^ t — 2/4 (a) and its inverse i 1— >■ a are of class C^. 

Thus in every Lorentzian frame we can represent our path in the form 

y = {r2{t),t) for all t G J, 

where t H> r2{t) is from some interval J into and represents the position of the 
mass as a function of time t in that Lorentzian frame. Clearly this representation 
is also of class and forms another equivalent parametric representation of the 
path but this representation is frame dependent. 

Most important parametrization of a path is with respect the proper time s 
of the moving mass mg. It is unique up to an additive constant and can be found 
from the formula 

(2.2) {dsf = {dy^f - {{dy^f + {dy^f + {dy^f) = {dtf - {dr^f. 
The condition 12.11 can be translated into 

dr2it) 



\r2{t)\ = 



< 1 = c for all t e J, 



dt 

that is velocity along any path of a point mass is less then the speed c of light. The 
proper time of a body carries over to general theory of relativity and thus it is also 
invariant under Lorentzian transformations. 

Maxwell established that waves in electromagnetic field propagate with velocity 
of light c. From considerations of Einstein and Rosen ^ follows that even distur- 
bances in gravity field should propagate with velocity of light. 

From results of Bogdan [1] and [2], Proposition 5.2, follows that if we consider 
the dynamics of n bodies interacting with each other by means of fields propagating 
with velocity of light, the equations of evolution are non-anticipating differential 
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equations and their solutions, not only depend on the initial conditions like in 
Newtonian mechanics, but also on the initial trajectory of the entire system. 

Assuming for instance that in a Lorentzian frame we are starting with n bodies 
whose initial trajectories 1 1— > yj{t), where j = 1, . . . , n, are known and we intend to 
observe the dynamics of evolution of the system for a period of time ti , and we can 
apriori estimate the bound v on velocities, and the bound A on the accelerations and 
the initial diameter 5 of the system, then the length of the interval of significance 
is at most, according to Proposition 5.2 of [5], 



Thus it is sufficient to know the initial trajectories of the system on the closed 
interval [a,0]. 

We should think about such trajectories as a postmortem record of the trajectory 
of some particular body from the system. It is clear that such trajectories would 
correspond to a time interval J that on the left is closed and on the right open or 
closed, finite or infinite. In any case it suffices to restrict ourselves to trajectories 
defined on intervals of the form J = [a, b) closed on the left and open on the right. 
The left end a of such time interval will be called a point of significance. Any 
time inside of the interval will be called a stopping time. 

Thus, if our trajectory t ^ r2{t) is of class C^, from continuity of the velocity 
w{t) = r2{t) and of acceleration w{t) on the closed interval [a,ti] follows that the 
following two functions 



are well defined for all stopping times ti € J, since the supremum of a continuous 
function on a closed bounded interval is attained at some point of that interval. 

For the sake of mathematical simplicity we shall consider only trajectories defined 
on the entire interval (— cx),cx)) = R. 

Definition 2.1 (Admissible trajectory). Assume that we are given a path of a 
point mass rriQ that in some Lorentzian frame has a representation in the form 
y — (r2{t),t), where the function r2{t) is from R into R^ and it has continuous 
derivatives up to order 3 and that for any stoping time ti J the kinetic energy 
and the acceleration r2{t) are bounded on the interval (— oo,ii). We shall say that 
such a function r2{t) represents an admissible trajectory. 

For the sake of logical completeness we should prove that every trajectory with 
a point of significance can be extended on the left to an admissible trajectory. Let 
us skip this for present and concentrate on admissible trajectories. 

Proposition 2.2 (Kinetic energy bound and velocity bound). Assume that a body 
having rest mass tuq moves along a trajectory r2 : i? — > R^. Let c denote the speed 
of light. For any nonnegative function k : R ^ R define function q : R ^ R by the 
formula 



Then for any t & R the following two conditions are equivalent 

• The kinetic energy of the body mo on the interval {—oo,t) is bounded by 



a = {S + 2vti)/{c~ v). 



(2.3) 



q{ti) — sup {1^(^)1 : u < ti,u e J} < c, 
A{ti) = sup{|?i;(u)| : u < e J} < oo. 




for all t e R. 



k{t). 
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• The velocity \v\ of the body mo on the interval (— oo,t) is hounded by cq{t). 



Proof. From Einstein's formula [8 , p. 22, the kinetic energy of mass mo moving 
with the velocity v is given by the formula 



moc^ ( — = — 1 



Thus the condition 

moc^ I ^ ^ 1 1 - ^(^) u<t,ueR 

\y/l-\v{u)\^/c^ J 

is equivalent to the condition 

\v{u)\ < cq{t) for all u < t, u £ R. 
This completes the proof. □ 



Notice that in the above proposition the quantity q{t) < 1 for all t £ R. 



Theorem 2.3 (Admissible trajectory is relativistic). The notion of an admissible 
trajectory does not depend on the Lorentzian frame. 

Proof. Assume that we have two Lorentzian frames S and S' . Assume that the 
frame S' moves away from frame S with constant velocity u. Assume that 1 1— >■ r2{t) 
represents an admissible trajectory in the frame S and a body with rest mass toq 
is moving along the trajectory. 

Without loss of generality we may assume that the frames S and S' are oriented 
so that the transformation of the coordinates y = (r, t) from S to S' is given by the 
formulas 

y'l = Vi 

2/2 = 2/2 

2/3=7 (2/3 - W2/4) 
2/4 = 7 (2/4 - U2/3) 

where 7 = (1 — v?)^^/"^ and and t/4 denote time in the respective frames. As 
before c = 1. 

First of all notice that the time interval (—00, 00) maps onto itself from frame S 
into 5". Indeed we have 

dy'i dy4-udy3 i i\ ^ n f u + 

- — — 7 = 7 (1 — uv) > 7 (1 — \u\) > tor all t = y^ £ R. 

dyi dyi 

Define function g by the formula 

9{t) = 2/4(2/4) for all t^y^eR. 
From Cauchy's mean value theorem we have 

g{t)^g{Q)=tg'{6)>t^ {l-\u\) for all i > 0. 
Thus y^ — g(t) 00 li t ^ 00 . Similarly 

c,(i) - c,(0) = tg'(6l) < <7 (1 - |w|) for all t<Q. 
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Thus 2/4 = g{t) — > — oo if f — > —oo. Since any continuous function maps an interval 
onto an interval the function g maps R onto R. 
Introduce a function / by the formula 

f{w) = (y=^ - l) for all w>0. 

Notice that the function / is nondecreasing and the kinetic energy of the mass mo 
moving along the trajectory can be represented as 

mofi\v\) 

where 

dys 

V = — = r2 

dj/4 

is the velocity of the body in the frame S. 

The velocity of the body in frame S' is given by 

I _ dy's _ dy3 - udy4 _ v-u 



dy'i dy4 -udys 1-uv' 

Thus we have the estimate 

I /, < 1^1 + 1^1 < + l^^l 
' l-\u\ - l-\u\ 

for all times in the initial interval {—oo,t). The quantity q{t) denotes the velocity 

bound on the initial interval. Thus the velocity is bounded on every initial 
interval (— oo,t') in the frame <S". Therefore its kinetic energy is bounded on every 
initial interval. 

Now let us consider the acceleration in the frame S' . It can be expressed as 

dv' v{l — V?) 
dy'i 7(1 - uv)^ 

in terms of quantities in frame S. Thus on every initial interval {—oo,t) we have 
the estimate 



dv' 



dy'i 



7(1 -H)3 

where A{t) is the bound on the acceleration in the initial time interval (—00, t) in 
the frame S. Hence the acceleration in the frame S' is bounded on every initial time 
interval (— oo,t'). 

Therefore the trajectory in the frame S' forms an admissible trajectory. □ 

3. Retarded time field 

Now consider any point (ri , in a fixed Lorcntzian frame and let (r2 (r) , t) 
denote a point on the path of the point mass with the property that light beam 
emitted from the trajectory will arrive at position n at time t. 

The time r is called the retcirded time. It must satisiy the relation 

|ri-r2(r)p-(i-r)2=0, 

which is preserved under Lorcntzian transformations. 

The following theorem establishes that the retarded time is well defined as a 
function of the variables (ri ,t) & R? x R. 
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Theorem 3.1 (The retarded tnne r is unique and forms a continuous function). 
Assume that we are given in a Lorentzian frame an admissible trajectory 1 1— > r2{t). 
Then for any point ri € and any time t £ R there exists a unique number t < t 
such that 

T = t-\ri~ r2{T)\. 

Moreover the map {ri,t) h- > r represents a locally Lipschitzian function on the 
space R^ x R. Thus T{ri,t) is continuous on R^ x R. 

Proof. For fixed ri G R^ and t £ R introduce a function / by the formula 

f{s) = t-\ri-r2{s)\ for ah s<t. 

The function / is well defined and maps the closed interval (—00, t) into itself. The 
function represents a contraction. Indeed 

\fis) - /(l)| = \{t-\r,~ r2is)\) -{t^ In - r2(l)|)| < {r^is) - r^^ 

v{x) dx\ < vi\s — s\ for all s,s<t, 

where vi = q(t) < c = 1 is the velocity bound corresponding to stopping time t. 
Therefore by Banach's contraction mapping theorem there exists one and only one 
solution of the equation t = /(t). 

Thus the map (ri , t) 1— > r is well defined in our Lorentzian frame for all points 
(ri,i) e R^ X R. 

To prove that the function r is locally Lipschitzian it suffices to prove that it is 
Lipschitzian on every open set of the form R^ x {—00, ti). To this end take any two 
points (ri, t) and (fi, t) from the domain of t such that t,t < ti. Let wi < 1 denote 
the velocity bound corresponding to our trajectory on the interval (—00,^1). 

To avoid unnecessarily complex notation denote by r and f the retarded times 
corresponding to the points (ri , t) and (fi , t) respectively. We have 

It - f I = |/(r) - f{f)\ = \{t-\r,~ r^ir)]) - {i - \h - r2(f)|)| 

< \t - i\ + In - h\ + \r2{T) - r2{f)\ 

= \t-i\ + \ri~h\ + \J^ r2{u)du\ 

<\t-i\ + \ri~h\+vi\T-f\. 

Taking the last term in the above inequality onto the left side and dividing by 
(1 — wi) both sides of the obtained inequality we get 

\T{n,t)-T{h,i)\ < + for all {ri,t),{h,i) & i?^x(-oo,ti) 

Thus the function t is continuous on the entire space R'^ x R. □ 

For a proof of Banach's contraction mapping theorem see, for instance, Loomis 
and Sternberg [13] page 229. 

Theorem 3.2 (An explicit formula for the retarded time function r). Assume that 
we are given in a Lorentzian frame an admissible trajectory t^r2 (t) . 

Take any stopping time ti £ R and let vi = q{ti) < c denote the corresponding 
velocity bound for t <ti. 
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Put So (ri ,t) ~ and define recursively the sequence 

SniriTt) — f{sn-i{ri,t)) for all n — 1,2,3,...; and rx ^ , t < t\, 

where f{s) =^ t — \ri — r2{t — s)\ for all s < t. 

The retarded function r is given by the formula 

T{ri,t) = lims„(ri,t) for all ri G and t £ R. 

n 

Moreover we have the following convenient estimate for the rate of convergence 

\T(ri,t) - sJri,t)\ < - — - — \t - \ri - r2(t)\\ for all ri e R^ and t < ti. 
1 - wi ' ' 

Proof. The proof follows from Theorem 4.7 page 37 of Bogdan or Theorem 9.1 
on page 229 of Loomis and Sternberg [13 . □ 



4. The fundamental fields associated with an admissible trajectory 

Now define the delay function T{ri,t) = t — T{ri,t) and notice that it satisfies 
the equation 

(4.1) r= |ri-r2(t-r)| for all t e R and ri e R^ . 

It is important to stress that in this paper we shall use the term field as syn- 
onymous with a function defined on a set of the space i?^. The space R"^ x R — R'^ 
is treated as a fixed Lorentzian frame. So we will use such expressions as scalar 
field, vector field, tensor field, etc., to describe the functions taking values in cor- 
responding spaces. Thus the functions r and T represent continuous scalar fields 
defined on the entire space i?**. 

Since the function T as difference of two continuous functions is continuous the 

set 

G = {{ri,t) eR^ X R: T{ri,t) > O} = T-^{Q,oo) 

as an inverse image of an open set by means of a continuous function is itself open. 
The set G consists of points that do not lie on the trajectory. 

By assumption the trajectory r2 has continuous derivatives f2{t) = w{t) and 
w{t), so we can define the vector fields 

ri2iri,t) = ri - r2{T{ri,t)), v{ri,t) ^ r2{T{ri,t)), a{ri,t) = w{T{ri,t)) 

for all {ri,t) € i?"*. Introduce the unit vector field e = ri2/T and fields u and z by 
the formulas ^ 

u = — and z ~ ■ on G. 

T {l-{e,v)) 

In the above (e, v) denotes the dot product of the vectors e and v. Since | (e, v) \ < 
jwl < c = 1 the vector field z is well defined on the set G. 

Definition 4.1 (Fundamental fields). Assume that we are given in a Lorentzian 
frame an admissible trajectory 1 1-^ r2 (t) . 

Define the time derivative w{t) — r2{t). The fields given by the formulas 

T, T, ri2 = ri — r2 o T, v = w o t, a — w o t for all (ri, t) Cz G 

and 

u—l/T, e — uri2, and z for all {ri,t) G G 
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will be called the fundamental fields associated with the trajectory r2{t). The 
operation o denotes here the composition of functions. 

The fundamental fields are continuous on their respective domains. This follows 
from the fact that composition of continuous functions yields a continuous function. 
Thus all of them, for sure, are continuous on the open set G of points that do not 
lie on the trajectory. 

Analogous fields defined by similar formulas on an open set G C i?^ appear in 
the problems involving plasma flows [3] and [1], or more generally flows of matter. 

We would like to stress here that the fundamental fields depend on the Lorentzian 
frame, in which we consider the trajectory. It is important to find expressions in- 
volving fundamental fields that yield fields invariant under Lorentzian transforma- 
tions. 

Lorentz and Einstein [8 , Part II, section 6, established that fields satisfying 
Maxwell equations are invariant under Lorentzian transformations. 

Our main goal is to prove that fields given by amended Feynman formulas and 
fields obtained from Lienard-Wiechert potentials satisfy Maxwell equations. We 
shall do this by showing that these fields are representable by means of fundamental 
fields and using the formulas for partial derivatives of the fundamental fields prove 
that such fields generate fields satisfying Maxwell equations. 

The following theorem represents the main pillar of the entire structure of the 
proof. Consider each of the formulas as bricks from which the pillar is constructed. 
If any one brick is rotten the whole structure of the proof will collapse. This is just 
a gentle warning to an impatient reader not to skip the computations involved. 

Introduce operators D = and Di = for z = 1, 2, 3 and V — {Di, D2, D^). 
Observe that 6i in the following formulas denotes the i-th unit vector of the 
standard base in that is 5i = (1,0,0), 62 = (0,1,0), ^3 = (0,0,1). 

Theorem 4.2 (Partial derivatives of fundamental fields). Assume that in some 
Lorentzian frame we are given an admissible trajectory t i~> r2{t). Define the time 
derivative w{t) — r2{t). For partial derivatives with respect to coordinates of the 



vector ri 


we have the 


following identities on the set G 


(4.2) 


D,T = 


zci where v — f2 t, 


(4.3) 


DiU = 


— ZU^Ci, 


(4.4) 


D,v = 


—CiZa where a — w t, 


(4.5) 


D,T = 




(4.6) 


Ae = 


—uzCiC + uSi + uzCiV where dt = (Sij), 


(4.7) 


D,z = 


-~z^ei{e, a) — uz^et -\- uz^a + uz'^Vi + uz^ei{v, v\ 


(4.8) 


VT = 


ze, 


(4.9) 


Vu = 


— ZM^e, 


(4.10) 


Vz = 


— z'^(e, a)e — uz^e + uz^e -f uz^v -f uz^{v, v)e. 
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and for the partial derivative with respect to time we have 



(4.11) DT = 1-z, 

(4.12) Du = zu^-u^, 

(4.13) Dt = z, 

(4.14) Dv — za where a = woT, 

(4.15) De — —ue + uze — uzv, 

(4.16) Dz = uz - 2uz^ + z^{e,a) + uz^ - uz^{v,v). 



Since the expression on the right side of each formula represents a continuous func- 
tion, the fundamental fields are at least of class on the set G. Moreover if the 
trajectory is of class C°° then also the fundamental fields are of class C°° on G. 



Proof. Proof of formula (|4.2p : Applying the operator Di to both sides of equation 
dm]) we get 

AT= (e,A(ri -ra)) = (e,(5y)) - (e,«(r))AT 
-e, - (e, v(T))(-Ar) = e, + {e,v)D,T 

yielding formula (|4.2p . 

Formulas (j4?3| . (|48l) . and (j4J| follow from formula (|42|) . 
Proof of formula (14. 4^ : 

DiV = aDiT — a{—DiT) ~ —Ciza 
Proof of formula (14. 5^ : 

D,T = D,{t - T) ^ - AT ^ -ze. 
Proof of formula (14. 6p : 

Ae = A[w(ri - r2)] = [Aw]ri2 + u[Diri - A?'2] 
= [Diu]u^^e + u[A?'i - A?'2] 
j^^^ = [— zu^ei]-u~"^e + u5i — u[At]w 

= —zuciC + U(5i — u[— zci]?; 
= —zucie + U(5i + uzciV 
= —uzae + + uzeiV 

Proof of formula (147)) : 

D,z = A(l - (e,«))-^ - {-l)z^{-D,{e,v)) 
= z'^[{e, Div) + (w, —uzeiC + + uzciv)] 
= z^[{e, —CiZa) + [v, ^zuCiC + u6i + uzCiv)] 

— -~z'^ei{e, a) — uz^ei{v, e) + uz'^Vi + uz^ei{v, v) 
= —z^ei{e, a) — uz^ei[l — z^^\ + uz^Vi + uz^ei{v, v) 

— ~z^ei{e, a) — uz^Ci + uz'^Ci + uz'^Vi + uz^ei{v, v) 
Formula (|4.10p follows from the formula (|4.7p . 



(4.18) 
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Proof of formula ()4.1ip : Applying the operator D to both sides of the equation 
(|4J|) we get 

DT = D\r^2\ = (e, ^['^12]) - -(e, Dr^) 

= -{e,v)DT = -{e,v){l - DT) = -{e,v} + {e,v}DT. 

The above yields 

— (e, v) 
1 - (e,u) 

Proof of formula (|4.12p : 

Du = DT-^ = {-1)T-^DT = (-l)u^(l - z) = zv? - . 
Proof of formula (|4TT3l) : 

Dt ^ Dif, ~T) ^ \ - DT = \ - {\ - z) = z. 
Proof of formula ()4.14p : 

Dv — aDr — za. 

Proof of formula (|4.15|) : 

De = D[uri2\ = [Du]ri2 + u[Dri2] 
= [zu^ — u'^\u~^e — u[Dr2\ 
= [zu — u\e — u[{Dt)v\ 
= zue — ue — uzv = —ue + uze — uzv. 

Proof of formula ()4.16p : 

Dz ^ D{1 - {e,v)y^ = (1 - {e,v)y^D{e,v) 
= z^{De,v) +z^{e,Dv) 

— z'^{[uze — ue — uzv], v) + z'^{e, za) 

— uz^{e, v) — uz'^{e, v) — uz^{v, v) + z'^(e, a) 

= - z^^] - uz^[l - z^^] - uz^{v, v) + z^{e, a) 

— uz^ — uz^ ~ uz^ + uz — uz^{v, v) + z'^{e, a) 
= uz^ — 2uz^ + uz ~ uz^{v, v) + z^{e, a) 

= uz — 2uz'^ + z^{e, a) + uz^ — uz^{v, v) 

□ 



5. Electromagnetic potentials satisfying wave and gauge equations 
induce field satisfying maxwell equations 

The following theorem represents the second main pillar of the structure. Each 
step one should prove in detail to see that no heuristic is introduced. The proof of 
this theorem one can trace back to the work of Lorentz [14] . 
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Theorem 5.1 (Wave and gauge imply Maxwell equations). Let on an open set 
G C i?^ he given two scalar fields (j) and S and two vector fields A and J. 

Assume that the fields (j) and A have second partial derivatives with respect to 
the coordinates of the point (ri ,t) G G and these derivatives are continuous on the 
entire set G. 

If these fields satisfy the following wave equations 
and Lorentz gauge formula 

ot 

on the set G, then the fields E and B defined by the formulas 

d 

E^-V(j)-—A and B = VxA for all {ri,t) e G 
will satisfy the following Maxwell equations 

(a) V-^ = 5, (b) VxE = -^B, 

(5.1) 

(c) V-B = 0, {d) VxB = —E+J 
on the entire set G. 



Proof. In this proof we will use the usual notation for dot product of two vectors 
F ■ H = {F,H) to be close to the customary notation used in textbooks. The 
following proof represents a reversed argument presented in Feynman-Leighton- 
Sands [TU], vol. 2, chapter 18, section 6. Namely we show that potentials, that 
satisfy wave and Lorentz gauge equations, induce a pair of fields E and B satisfying 
Maxwell equations. 

Since for any vector field F, which is twice continuously differentiable, we have 
V • (V X F) = 0, the vector field B must satisfy Maxwell's equation (c). 

Since the differential operators V and D commute, due to the assumption that 
the fields (p and A are of class , and since from the Lorentz gauge formula follows 
that W ■ A = —D(j>, we must have 

V • = V • (-V?!) - DA) = ~D\/ ■ A = -V^(/. + D'^(t) = S. 

Therefore the equation (a) is satisfied. 

Since V x (V/i) = for any scalar function, that is twice continuously differen- 
tiable on the set G, and for any such vector function F defined on the open set G 
we have V x DF = DV x F, we get 

V X S = V X - DA) = -DV X A = -DB 

and this means that Maxwell's equation (6) is satisfied. 

Finally for the last equation. Maxwell's equation (d), from the algebraic identity 

V X (V X A) = V(V • A) - V^^ 

and Lorentz gauge equation and the wave equation 

y ■ A + D4) = Q, V'^A-D'^A^-J 
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we get 



V X B-DE = \I X (V X ^) - D{-\I<f)- DA) 
= V(V • A) - V^A + DV(t) + D'^A 
= V[(V • A) + D(t)\ - {V^A - D^A) = J. 



□ 



Notice that in the above theorem the field S corresponds to the field g of charge 

density and J to the field j of intensity of charge flow. Since we redefined the 
electrostatic constant eo to satisfy the condition Aneo = 1 our new unit of charge 
became equal to Aw [coulomb] . 

It is worthwhile to mention also the following proposition. 



Proposition 5.2 (Equation of continuity for flow of charge). If the electromagnetic 
potentials (j) and A have continuous partial derivatives of order 3 on an open set 
G C i?**, then from the Lorentz gauge equation 

D(I> + V-A = Q for all {ri,t) G G 

follows the equation 

DS + W-J = for all {ri,t) G G 
of continuity for charge flow. 

Proof. Indeed introduce the D'Alembertian operator 



Prom our assumptions follows the commutativity of the differential operators 



D, V, 



□2 



when acted on flelds (p and A, respectively. Hence we have 



a^D = Da^ and a^\/ = \/a^. 



The above implies 



DS + V ■ J = D tf^ + V • D'^A = n^{D(t) + V ■ A) =0 



and this is the equation of continuity of flow of charge. 



□ 
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6. LlENARD-WlECHERT POTENTIALS AND WAVE AND GAUGE EQUATIONS 

In this section we will deal only with Lienard-Wiechert potentials 

(/) = uz and A = uzv. 

The author in the proof of the following theorem used DOS based interactive 
computer software called DERIVE allowing to input user data files and perform 
on them interactive operations. More derails about DERIVE and actual programs 
used in computation we shall present in the last section of the paper. 

Any program, that can symbolically manipulate polynomials of an arbitrary 
number of variables, can be used to verify the proof. For a reader, who would like 
personally to verify the proof by hand or another software system, we present here 
just the key intermediate formulas as computational mile stones. 

In the following formulas the vector field d denotes the composition of the func- 
tion ^r^it) with the retarded time field r. 

Theorem 6.1 (Homogenous wave and gauge equations are satisfied). For any 

admissible trajectory r2 {t) of class the two homogenous wave equations 

{W^ - D^)[uz]=0, {V^ -D^)[uzv]=0, 
and Lorentz gauge equation 

V • [uzv] + D[uz] = 
are satisfied on the entire set G = {{ri,t) G : T{ri,t) > 0} . 

Proof. We have the following identities 

V^M = D^[uz] = + uz'^ie, a) + u^z^{e. a) + Suz'^ie. a)^ + 3u^z^ 

+ 3u^z^ + 3u^z^{v,vf +6u'^z''{e,a) +6u^z^{v,v) 

— 3v^z^{v, a) — Qu^z'^{e, a) — &u^z^{v, v){e, a) — Qu^z'^ 

— &v!^z^{v,v) — u'^z^{v,v) 
Thus the first wave equation takes the form 

V^[uz]- D^[uz] = 0. 

Similarly for the second wave equation we have the formulas 

V'^luzv] = D'^luzv] = 

= +uz^d + uz^{e, a)v + u^z^a + u^z^{e, a)v 

+ 3uz^{e, a)a + 3uz^{e, a)'^v + Su'^z'^a + Su'^z^v 

+ 3u^z^{v,v)'^v + 3u^z^v + 6u'^z'^{e,a)v + 6u^z'^{v,v)v 

— u'^z^{v, v)v — 3m^2;^(w, a)v — Su^z^lv, v)a — A^u^z^a 

— &u^z'^{e, a)v — 6u'^z^{v, v){e, a)v — du^z^v — 6u^z^{v, v)v. 
Computing terms of Lorentz gauge equation we get 

V • [uzv] = (V, uzv) = +u'^z^{v, v) — u^z^ + v?7? — uz^{e, a) 
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and 

D[uz] = ~u'^z^{v,v) + u^z^ - u^z^ + uz^{e,a). 

Hence 

{V,uzv) + D[uz] = 0. 

□ 

Lienard-Wiechert potentials generate a solution to a homogenous system of 
Maxwell equations. Indeed from previous theorems we get the following corollary. 

Corollary 6.2 (Lienard-Wiechert potentials and Maxwell equations). Let r2{t) he 
any admissible trajectory and let G denote the open set of points that do not lie on 
the trajectory, that is 

G = {(ri,t) e R'^ : T{ri,t) > O} . 
Define the fields by formulas 

(j) — uz, A — uzv, S ^ 0, J = for all {ri,t) G G. 
Then the fields E and B defined by the equalities 

d 

E ^ -V<j) - —A and B^VxA for all {ri,t) e G 

will satisfy the following homogenous Maxwell equations 

d d 
V-E^O, y X E+—B ^0, V-B = 0, yxB-—E = 
ot at 

on the set G. 

7. Proof of the amended Feynman's formulas 

As before let r2{t) be any admissible trajectory of class G^ in some Lorentzian 
frame and let 

G = {{ri,t) € : r(ri,t)>0} 
be the set of points not lying on the trajectory. Let cf) = uz and A = uzv denote 
the Lienard-Wiechert potentials corresponding to the trajectory r2. 

We remind the reader that all intensities of the considered fields are per unit of 
charge, thus if we are dealing with a charge q the corresponding potentials will take 
the form (j) — quz and A — quzv and similar adjustments should be made with the 
fields E, B. 

Feynman using his formula for the electromagnetic potentials, see vol. 2, Feynman- 
Leighton-Sands [T^, chapter 15, page 15.15, had derived on page 21.9, of chapter 
21, a formula for electromagnetic potential of a moving point charge. A formula, 
which did not agree with the Lienard-Wiechert formulas. 

Nevertheless on page 21.11 in the last sentence he suggested that perhaps the 
reader may take his word for it, that Lienard-Wiechert potentials yield the same 
field as his formulas in the case of a moving point charge. In the footnote he 
mentioned that he has done that, but it took him a lot of time and paper. 
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This little footnote provided the author with the main impetus to try to verify the 
computations. Feynman did not leave any hints concerning possible simplifications 
of the computations or that quantities involved should be considered as fields. 

The author taking Feynman's suggestion performed by hand, step by step, com- 
putations presented in the following proof. These computations have provided an 
essential hint how to develop and to structure the proofs of the theorems presented 
in this paper. 

Thus the credit has to be given with gratitude to Feynman for discovery of this 
gem stone of Einstein's Special Theory of Relativity and consequently of Maxwell's 
Theory of Electrodynamics. 

To verify in a more efficient way that the electric fields E and Ef coincide, we 
provide a computer program in the last section of this paper. However to show 
that the magnetic fields B and Bj are equal a direct computation seems to be more 
efficient. 

To Feynman are due immense thanks and gratitude beyond words! Here are the 
computations as they might have looked under Feynman's hand. 

Theorem 7.1 (Equality of amended Feynman and Lienard-Wiechert fields). As- 
sume that in some Lorentzian frame we are given an admissible trajectory r2{t) of 
class . 
Let 

G ^ {iri,t) e : T(ri,i)>0} 

denote the set of points not lying on the trajectory. 

Let u and e denote the fundamental fields J[ ) generated by the trajectory 7*2 (i). 
Define fields 

E = -Wet)- DA and B = W x A, 

Ef ^ u^e + u^^D{u^e) + D'^e and Bf^exEj, 

Then we have the equalities E — Ef and B ~ Bf on the set G. 



Proof. From the formulas (I4.10p and (|4.9p we get 

V(t> = V(zu) = uVz + zVu 

= u[—z^{e, a)e — uz^e + uz^e + uz^v + uz^{v, v)e] + z\—z'u^e\ 

r 3/ \ 23|22|22|23/ \1 22 

= [— uz (e, a)e — u z e + u z e + u z v + u z [v, v)e\ — u z e 
— —uz^{e, a)e + u^z^v — u^z^e + v?z^{v, v)e 

and 

D[uv] = [Du]v + uDv — [zu^ — u'^]v + uza 

2 2 

= u zv — u V + uza 
— uza — u^v + u^zv 
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DA = D[zuv] = [Dz]uv + zD[uv] 
= [Dz]uv + zD[uv] 

= [uz — 2uz^ + 2;^(e, a) + uz^ — uz^{v, v)]uv 

+ z[uza — u^v + u^zv] 
= u zv — 2u z V + uz {e, a)v + u z v 

— v?z^{v,v)v + uz'^a — v^zv + u^z^v 
= uz^a — v?z^v + uz^{e, a)v + u^z^v — u'^z^{v, v)v 

thus 

E = -Vcf) - DA 

= — [—uz^{e, a)e + i^z^v — u^z^e + u'^z^{v, v)e\ 

— [uz^a — v^tP'v + U7^{e, a)v + u^z^v — u'^z^{v, v)v] 
= — [uz^a — uz^{e, a)e + uz^{e, a)v 
— v?z^e + v?z^{v^ v)e + u^z^v — u'^z^{v, v)v] 

Now consider the second term of the expression Ef. We have 

u-'^D[u^e\ = u-'^[Du^]e + uDe 

= u~^2u[Du\e + u[—ue + uze — uzv] 
= 2[u'^z - v?\e - M^e + ?i^ze - -u^ zv 
= 2v?ze — 2v?e — ii^e + u^ze — u^zv 
= Su^ze — 3u^e — u^zv 
= — 3u^e + Su^ze — u^zv 
= 2u^{z — l)e + u[uze — ue — uzv] 
= 2u'^ze — 2u^e + u^ze — u^e — u^zv 
= Su'^ze — 3u^e — u^zv 

For the third term of the expression Ef we get 

D^e = D[De] = D[—ue + uze — uzv] = D{u[z{e — v) — e]} 
= {Du}[z{e — f ) — e] + uD[z{e — v) — e] 
= {x?z — u'^}[ze — zv — e] 

+ u[{Dz){e -v) + zD{e - v) - De] 
= 2u^e — Su^ze + u^zv — uz^a + uz^{e, a)e — uz^{e, a)v 

+ iv^z^e — u'^z^{v, v)e — u^z^v + u'^z^{v, v)v 
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Compute the field Ef 

Ef = it^e + u-^D{u'^e) + D'^e 

— u^e + [— 3u^e + Su^ze — u^zv] 

+ [2u^e — Zu^ze + v? zv — uz'^a + uz^{e, a)e — uz^ {e, a)v 

+ z^e — u^z^ {v, v)e — z^v + z^{v, v)v 
= —uz'^a + uz^{e, a)e — uz^{e, a)v 

+ u^z^e — V? {v^ v)e — u^z^v + u'^z^{v, v)v 
and compare with 

E = —[uz^a — uz^{e, a)e + uz^{e, a)v 

— u^z^e + u'^z^{v, v)e + u^z^v — u^z^{v, v)v]. 
Now the magnetic field from Feynman's formula is 

Bf = e X Ef —ex [—uz'^a — uz^{e, a)v — u^z^v + u'^z^{v, v)v] 
— —uz'^ e X a — uz^{e, a) e x v — u^z^ e x v + v? z^{v, v) e x v 

and since 

V[u2] — —uz^{e, a)e + u^z'^v — u^z^e + u'^z^{v, v)e, 
\7 X V ~ ~-z e X a, 
we have from Lienard-Wiechert potentials that 

B = VxA = Vx [uzv] = [uz] V X V + V[uz] x v 
= [uz]{—z e X a) + {—uz^{e, a)e + z^v 

— u^z^e + u^z^{v, v)e) x v 
= —uz^ e X a — uz^{e, a) e x v — u^z^ e x v + u'^z^{v, v) e x v 
Clearly we have the equality for both fields 

E^Ef and B = Bf. 
Eureka! A salute is due to Feynman! □ 

Now let us collect the formulas representing the fields E and B by means of the 
fundamental fields (|4.ip into one convenient statement. These formulas can be very 
useful for computations especially with the help of a digital computer. 

Notice that all the fundamental fields generated by the admissible trajectory are 
defined by means of algebraic operations and operation of composition of functions, 
with the exception of the retarded time field r, but Banach's fixed point theorem 
(j3.2p provides a fast converging algorithm for computing this function. 

Proposition 7.2 (Representation of amended Feynman's fields by fundamental 
fields). Assume that in some Lorentzian frame we are given an admissible trajectory 
r2(t) of class of a moving point mass rriQ. Let 

G = {iri,t) e : T{ri,t) > 0} 

denote the set of points not lying on the trajectory. Define fields 

E = u^e + u-^D{u^e) + D'^e and B^exE, 
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Then the fields E and B have the following representation in terms of the funda- 
mental fields (^TTp generated by the trajectory r2{t). 

E = —uz^a + uz'^{e, a)e — uz'^{e, a)v 

23 23/ \ 23 23/ \ 

+ u z e~u z {v,v)e — u z v + u z {v,v)v, 
B = —uz^ e X a — uz^{e, a) ex v — v^z^ e x w + u'^z^{v, v) e x v, 

8. Modified theory of generalized functions 

Now let us prove that any pair {E, B) of fields of class on an open set G 
in R!^ that satisfies the homogenous system of Maxwell equations must satisfy the 
homogenous wave equations on G. 

It will be convenient here to use the theory of generalized functions, called also 
distributions, originally introduced heuristically by Dirac in his works on Quantum 
Mechanics and put on precise mathematical footing by L. Schwartz |15j and Gelfand 
and Shilov JJ^ . We do not introduce here any notions from the theory of Lebesgue- 
Bochner integral. We base the proofs only on elementary properties of continuous 
and differentiable functions that can be found in any book on Calculus. 

We will need to make a few modifications for the sake of our problem. We need 
distributions defined over any open set G C i?^ and having their values in the vector 
space R'^. 

Let Vk denote the space C^{G, R'^) of infinitely differentiable functions having 
compact supports contained in G and values in R'^. 

On the space Vk we introduce the usual sequential topology. A sequence of 
functions gn € T)/- converges to a function g g 2?^ if it converges uniformly together 
with all its partial derivatives D'^g, on every compact subset K of G, to the function 
g, that is for every a the following sequence converges uniformly on K 

D'^g^ix) ^ D"g{x) 

where a = (fci, • • • , k^) denotes the multi-index of a partial derivative 

L»" = D^i . . . 

with kj 0, 1, 2, . . . and Dj = g|- where j = 1, 2, 3, 4. 

Definition 8.1 (Generalized vector valued functions). Let 2?[. denote the space of 
all linear continuous real functionals on Vk- Convergence in Vf. will be understood 
as pointwise convergence. The space D'f, will be called the space of generalized 
vector valued functions or the space of vector distributions. 

We shall use the following equivalent notation for such a functional 
(8.1) /(<?) = (/,<?) = / f{x) ■ g{x) dx for all g G Vu. 

JG 

In the above y ■ x denotes the dot product, also called the scalar product, of two 
vectors y,x S i?*^. 

Any continuous function / on the set G generates by means of the above integral 
formula a vector distribution. 
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Notice that the space I?^ is linearly and topologically isomorphic with the Carte- 
sian product of k copies of the space 

V[x ■■■ xV[ = 

Indeed, if for every m = 1, . . . , fc we define maps Pm : -R — >■ -R'^ by the condition 
Pm{t) = a; {xm = t and Xj = if j ^ m}, 
where x = {xi, . . . ,Xk) € R'^, then the functional 

fmig) = / fix) ■ Pmigix)) dx for all gGVi 

Jg 

are well defined and represent elements of the space 2?^. Thus the element 

(/!,••• ,fk)€{V[)''. 
Conversely define maps : i?'^ — >■ for m = 1, . . . , by the formula 
P^{x) = Xm for all x G R!^, m= 1,. . . ,k. 

Then the formula 

fig) = E / fmPMx)) dx for ah geVk 

yields a linear continuous functional on the space Vk. It is easy to verify that the 
transformation Q : ^ (D^)*' defined by 

(/i, •••,//=) 

is indeed a linear and topological isomorphism of the two spaces. 

Now notice the following fact. 

Proposition 8.2 (Imbedding of continuous functions into V/. is one-to-one). Given 
two continuous functions fi and /2 on the open set G such that they generate the 
same vector distribution, that is 

/ fi{x)g{x)dx= / f2ix)gix)dx for all g€Vk, 
Jg Jg 

then they coincide 

fi{x) = f2ix) for all xgG. 



Proof. Indeed, from linearity of the map / (/, g) and the previous isomorphism 
it is sufficient to prove that for any real continuous function / such that 

/ f{x)g{x) dx = Q for all g 
Jg 

follows that / = on G. 

Assuming that this is not true then at some point a;o € G we have /(xo) ^ 0. 

We may assume without loss of generality that 25 = /(xo) > otherwise we would 
consider the function — /. From continuity of / follows that there is a rectangular 
neighborhood V C G oi xq such that 



fix) > S for all xgV. 
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There exists a nonnegative function g of class C°° with support in the set V 
with integral jQg{x) dx = 1. Thus for such a function we would get 

/ f{x)g{x)dx= i f{x)g{x)dx> / 5g{x)dx = S>0. 
Jg Jv Jv 

A contradiction. So all we have to do is to show that there exists a function g 

having the above properties. To this end consider a nonnegative function go{t) 

defined by the formula 

go{t) = ae-^/(^-*') if \t\ < 1; and go{t) = if \t\ > 1, 

where the constant is selected so that j\go{t)dt = 1. It follows from the above 
formula that the function go is infinitely differentiable at every point except perhaps 
atf = +lort=— 1. One can prove that at these points the one-sided derivatives 
exist and they arc equal. Thus the function go is of class C°°. 

Now consider the function 

giix) = goixi)go{x3)go{x2)go{x4) for all x = {xi,x2,X3,X4:) G 

with the integral, over the cube in representing its support, equal to 1. For the 
fixed xq G G and sufficiently large n we see that the function given by the formula 

gn{x) = n^g{n{x — xq)) for all x G 

will satisfy our requirements. □ 

Any linear continuous operator H : Vk Vf. generates a linear continuous dual 
operator H' : — )■ 2?^ by the formula 

{H'f,g) = {f,Hg) for ah geV^. 

The dual operator corresponding to scalar multiplication g i-^ Xg is scalar mul- 
tiplication / I-)- A / as follows from the above definition. 

Definition 8.3 (Generalized differential operator). By a generalized partial 

derivative = acting onto the m-th component of f G 2?^ 

Di^mifl, ■ ■ ■ , fk) = • • • 1 Di fm, • • • , fk) 

we shall understand the dual operator to the operator acting onto the m-th compo- 
nent of g GVk by the formula 

(-l)A,m(5'l>---,fl'fc) = • • • ; (-l)A5m, • • • iSfc)- 

In the case when the set G represents a Cartesian product of open bounded 
intervals and the vector function / is continuous together with Difm, the above 
formula can be easily verified through iterated integral and integration by parts. 
In this case the ordinary partial derivative will produce the derivative in the sense 
of distributions. So it is natural to extend this property to vector distributions. 

Definition 8.4 (Weak and strong partial derivatives). Now every distribution and 
in particular every continuous function is differentiable in the sense of distributions. 
If it happens that the distributional partial derivative is representable by means of 
a continuous function, such a function is called a weak derivative. // a vector 
function has a continuous partial derivative such a derivative is called a strong 
derivative. 
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It is not obvious that weak and strong derivatives coincide in the general case of 
an arbitrary open set G in i?"*. This calls for the following theorem. 



Theorem 8.5 (Weak and strong partial derivatives coincide). Assume that f = 
(/i, . . . , fk) is a vector-valued function on an open set G C R'^ and that its m-th 
component has a continuous partial derivative Dif„i. 

Then this derivative coincides with the derivative in the sense of distribution, 
that is the weak derivative coincides with the strong one. 



Proof. To prove this fact for general open sets in notice that for every point 
X & G there exists a neighborhood V{x) C G in the form of the Cartesian product 
of open intervals. Restricting our test functions to functions with support in V{x) 
will yield that on such a domain the weak and strong partial derivatives coincide. 

Since every open set in is a union of a countable number of compact sets, 
from the collection V{x) {x E G) one can extract a locally finite countable cover of 
the set G. 

Now using partition of unity theorem, (for reference concerning this theorem see 
for instance Gelfand and Shilov [11], vol. 1, Appendix to Chapter 1, Section 2,) we 
can prove this theorem for any open set G C R^. □ 



Now let A denote the algebra of all linear transformations H : PJ. — > T)'^. The 
multiplication of two transformations H and S in the algebra A is defined as com- 
position H o S. 

Since any two partial derivatives commute and they commute with scalar mul- 
tiplication operators, the above definition 18.31 permits one to represent any polyno- 
mial of the operators Di^„i as a dual of some operator acting in the space Pfc. It 
follows from the above definition that any operator H, generated by a homogenous 
quadratic polynomial p , is formally self dual 

{Hf,g) = {f,Hg) for all geVk- 

Now let us consider a particular space of interest £ —T)'^ x Pg. 
Let Aq denote the algebra of all linear transformations from into itself with 
multiplication defined again as a composition of transformations. 

Proposition 8.6 (Matrix representation of operators on £). A transformation H 
belongs to the algebra A if and only if there exists a 2 x 2 matrix such that 

"f'l ^ ^' 
72. 

where transformations Hij belong to algebra Aq. This correspondence establishes 
isomorphism between the algebra A and the algebra of all 2 x 2 matrices with ele- 
ments in the algebra Ao- 



(8.2) 



H{f) 



Hii,Hi2 




7i" 


H21,H22 




A 



for all f 



Proof. The proof is evident and we leave it to the reader. 



□ 



A transformation H with matrix of the form aiji, where a,;j G R and / is the 
identity map in Aq, will be called scalar transformation. A scalar transformation 
H can be simply represented as 



H{f)^ 



ail 


ai2 




7i" 


021 


a22 




A 



for all f e£. 
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Notice that the identity transformation, the unit element / in the algebra A, can 
be represented as a scalar transformation with matrix 

1 0' 
1 

Of a particular interest are the following two scalar transformations that are 
inverse to each other 



J 



-1 




and J 



1 

-1 



J. 



A transformation H E A is called simple if its matrix representation is of the 
form 

'h 0" 



(8.3) 



h 



where h e Ao- 



We shall use a shorthand notation when dealing with an action by a simple 
transformation. We shall write h{f) or just hf instead of 



(8.4) 



h 
h 



for all 



= .fe£ 



Proposition 8.7 (Commutativity of J with any simple transformation). Every 
simple transformation H d A commutes with transformation J that is 

HoJ = JoH 

In shorthand notation we can write hj — Jh if no confusion is possible. 

Proof. The proof is obvious and we leave it to the reader. □ 

Compare the above constructions with the arguments presented in Kaiser |12) . 
page 207. 

Theorem 8.8 (Homogenous Maxwell equations imply homogenous wave equation). 
Let E and H be two vector distributions from the space Pg. // they satisfy the 
homogenous .system of Maxwell equations 

V X E = -DB, V ■E = 0, 

VxB = +DE, V-B = 0, 

then they satisfy the wave equations 

(V^ - D^)E = 0, 
(V^ - D^)B = 0. 



Proof. Let 

■' [B 

Then Maxwell's equations can be written as 

(a) {\/x)f = DJf, (&) (V-)/-0. 
Applying operator J^^D — — JD to both sides of equation (a) we get 
(8.5) -JD{Vx)f^D^f. 
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On the other hand applymg the simple operator (Vx) to both sides of equation 
(a) and using the commutativity of the operators J, £>, (Vx) we get 

(8.6) (Vx)2/-(Vx)Ji?/ = Ji?(Vx)/. 
Since it is easy to verify that 

(8.7) (Vx)2/ = V((V-)/)-VV for all feS, 

observing that the terms on the right side represent differential operators generated 
by quadratic homogenous polynomials and the above identity holds for any test 
function from V^, by Maxwell's equation (5) we get that the solution / of the 
system of Maxwell equations must satisfy the equation 

(Vx)^/ = -VV. 

Thus from equations 18. 5[ 18.61 and 18.71 follows that the solution / of the system of 
Maxwell equations must satisfy the wave equation 

and this completes the proof. □ 



Proposition 8.9. Assume that E and B are two, twice continuously differentiable, 
vector fields from an open set G in the space into B? satisfying the homogenous 
system of Maxwell's equations 

VxE = -DB, V-E = 0, 

V X B = +DE, V • S = 0, 

Then these fields satisfy the homogenous wave equations 

V'^E- D'^E = 0, V^B - D'^B = 0. 



Proof. The proof follows from the previous theorem and from the fact that for 
continuous functions with continuous partials their partial derivatives in the sense 
of generalized functions coincide with the corresponding ordinary partial derivatives 
that is the weak derivatives coincides with the strong derivatives, see Theorem 
(18311 . □ 

9. Maxwell's mathematical experiment 

At the beginning of the 19-th century the magnetic and electric phenomena 
were considered as independent. Ampere's law that electric current can produce a 
magnetic field was the first breakthrough in this field. 

Around 1830 M. Faraday confirmed validity of Ampere's law. Further he dis- 
covered that changing magnetic flux can induce an electric current. He established 
the principle of electric induction leading to transformers of alternating current. 

In 1846 Faraday wrote a paper Thoughts on ray vibrations. In his imagination he 
saw interaction of gravitational, electric, and magnetic forces along the force lines, 
which connected particles and masses together, contrary to prevailing Newtonian 
view of actions at a distance. Faraday in his work suggested that light itself is some 
kind of vibration in the lines of force. He postulated that vibrating charge may 
produce vibrations in the lines of force. 
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J. C. Maxwell was inspired by Faraday's work and after augmenting Ampere's 
law he discovered the specific velocity c associated with electromagnetism that is 
now known as the speed of light and proved that electromagnetic waves in free 
space propagate with velocity c. 

We can perform this experiment now in a more general setting and a greater 
generality. Maxwell had done it for harmonic waves in electromagnetic field. 

Let us assume that we are working in some Lorentzian frame and we have some 
physical quantity u, it could be a component of the electric or magnetic field, gravity 
field, or any other quantity, that satisfies the homogenous wave equation in some 
open set G G R*. We can select a rectangular subset of G, so let us assume that 
G is in rectangular form, that is, it is representable as a Cartesian product of 4 
intervals. The wave equation in standard g-m-s units for quantity u will look as 
follows 

(9.1) V^u-— — w = for all (r,t)eG. 

We may assume without loss of generality that the origin of our Lorentzian frame 
is in the center of the 4-dimensional cube G. Take any function f{t) of class G^ on 
R and consider a wave that this function generates by means of the formula 

u{yi,y2,y3,t) ^ f{yi - vt) for all {yi,y2,y3,t) E G. 

At time t — the shape of the wave has graph represented by function the t !—> f{t). 
To visualize this wave we may assume that the wave looks like tsunami wave with 
the crest at i = 0. At any time t the crest is at position yi = vt. Thus the wave 
moves with velocity v. 

We want to find when such a wave will satisfy our wave equation. Notice that u 
does not depend on 2/2 and y^ thus 

V^u = f"{y,-vt) for all (yi, 2/2, ys, G G. 

and 

—u = v^f"{yi~vt) for all (yi, y2, 2/3, i) e G. 

The wave equation (|9.ip will be satisfied for any function of class G^ if and only if 
= v'^ that is w = +c or v — — c. Thus such waves propagate with velocity c. This 
is the essence of Maxwell's mathematical experiment. Hence we have the following 
corollary 

Corollary 9.1 (Maxwell's experiment). Assume that in a Lorentzian frame we 
have some physical quantity u satisfying the homogenous wave equation \9.1\) . then 
waves of that quantity propagate in free space with velocity of light c. 



10. Einstein's formulas for transformation 
of electromagnetic field 

Lorentz [14^ and Einstein [S] , Part 2, section 6, established that a pair of fields 
E^{Ei,E2,E3) and B={Bi,B2,B3) 
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that satisfy homogenous system of Maxwell equations in an open set G transforms 
as a part of an antisymmetric tensor of second rank in Lorentzian space-time x R. 
The matrix of this tensor looks as follows 

+Ei +E2 +E3 

— El — B2 

—E2 —B^ +Bi 

—£"3 +-B2 —Bi 

For formulas for transformation of electromagnetic field under Lorentzian change 
of coordinates see Feynman-Leighton-Sands [10], vol. 2, chapter 26. In particular 
formulas for Lorentz transformations of the pair of fields E and B on page 26-9. 
Compare also with Dirac [7], section 23, page 42. 

Following [TU] assume that our orthogonal frames of reference are S and S' and 
the frame S' moves along the first coordinate axis of S with velocity w — [wi, 0, 0). 
Then the transformation of the components of the vector fields 

E = {Ei,E2,E3) and B = [Bi, B2, B:i) 

will look as follows 

E[ ~ El B[ ^ Bi 

E'2 = "/{E + WX B)2 B'2= j{B - ?« X E)2 

E'^ = "f{E + WX B)3 = 7(B -wx E)^ 

where 7=1/ \/l — 
Since 

w X B = w X [e X E) ~ {w ■ E) e — [w ■ e) E 
we get final formula involving just the vector fields E and e 

E'l^Ei 

E2 = 7(-E'2 + [w ■ E)e2- {w ■ e) E2) = l{E2 + wiEi 62 - wiCi E2) 
E3 = liEs + iw ■ E)e3~ {w ■ e) £'3) = 7(£;3 + wiEi 63 - wici E3) 

11. Bogdan-Feynman Theorem for a moving point mass 

From the theorems of the previous sections we can conclude that the following 
theorem is true. Please notice that the physical nature of the fields is completely 
irrelevant. 

The important fact is that we are working in a Lorentzian frame, in which we are 
given an admissible trajectory, see definition (|2.ip . This trajectory alone generates 
in a unique way the system of the fundamental fields by means of which we are able 
to construct fields that are preserved by Lorentzian transformations. 

As we have proved earlier the notion of an admissible trajectory is invariant 
under any Lorentzian transformation from one frame to any other frame moving 
with a constant velocity. This fact is one of the consequences of Einstein's formula 
on kinetic energy of a point mass. 

We remind the reader that we are working in a fixed Lorentzian frame. The 
position vector is denoted by ri S R^ and time hy t G R. The position axes are 
oriented so as to form right screw orientation. The units are selected so that the 
speed of light c = 1. 
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Partial derivatives with respect to coordinates of ri are denoted by Di, D2, -D3 
and with respect to time just by D. The gradient differential operator is denoted 
by V = {Di,D2, D3) and the D'Alembertian operator by = - D^. 

Theorem 11.1 (Bogdan-Feynman Theorem). Assume that in a given Lorentzian 
frame the map t r2 (t) from R to represents an admissible trajectory of class 

. Assume that G denotes the open set of points that do not lie on the trajectory. 
All the following field equations are satisfied on the entire set G. 

Consider the pair of fields E and B over the set G given by the formulas 

E = u^e + u^^D{u'^e) + D^e and B = exE 

where u and e represent fundamental fields {-j.l^ associated with the trajectory r2 (t) . 

Then this pair of fields will satisfy the following homogenous system of Maxwell 
equations 

VxE=-DB, V-E = 0, 

VxB = +DE, V-B = 0, 

and the homogenous wave equations 

D^E = 0, D^B = 0. 

Moreover Lienard-Wiechert potentials, expressed in terms of the fundamental 
fields as A — uzv and (p — uz, satisfy the homogenous system of wave equations 
with Lorentz gauge formula 

□^^ = 0, tf0 = O, V ■ A + D4)^Q 

and generate the fields E and B by the formulas 

E = -V0 -DA and B = W x A, 

Finally we have the following explicit formula for the field E in terms of the 
fundamental fields 

E = —uz'^a + uz^{e, a)e — uz^{e, a)v 

23 23/ \ 23 23/ \ 

+ u z e — u z {v,v)e~u z v + u z {v,v)v. 

As a consequence of the above theorem the components of the quantities E, B, 
A, and (j) propagate in the Lorentzian frame with velocity of light c = 1. 

12. Einstein's illusive gravity field 

Matter in the universe is distributed in chunks like atoms, particles, molecules, 
or stellar objects, like planets, stars, galaxies, etc. So it is natural to consider a 
system on n bodies carrying charges qj and having rest masses moj ■ 

Einstein using general theory of relativity proved that waves in gravity fields 
should also propagate in space with velocity c like electromagnetic waves. For 
reference see Einstein and Rosen [5]. We shall assume that the field E given by 
Newton-Feynman formula has this property. We want to find out how to develop 
dynamics in such fields. 

Let us assume that we are working in a Lorentzian frame with units selected 
so that the speed of light c = 1 the electrostatic constant in free space satisfies 
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the condition 47reo = 1 and the unit of mass is selected so that the gravitational 
constant G = 1. 

This means that if we keep meter as our unit of length the unit of time is 
approximately 3.3 nanoseconds and the unit of mass is about 3.871 metric tons and 
the unit of charge is 47r of coulombs. 

Let us assume that we have a system consisting of n bodies indexed by j, k = 
1, 2, . . . , n with rest masses moj and charges qj. We want to find the equations of 
motion for a time period starting at time t = s and ending at time t — e. Assume 
that the bodies during this time period [s, e] — I move without collision. 

We shall follow here the constructions developed in Bogdan [1: . Now consider the 
trajectories. Assume that we know the formulas for the trajectories of the bodies 
up to time t — s. We can setup the equations of motion of these bodies in their 
own force fields. To this end let rj{t) be the trajectory of the j-th body. Let tjk{t) 
be the retarded time for a wave travelling at the speed of light to reach from the 
k-th trajectory the j-th trajectory at the point rj(t) at time t. 

This function must satisfy the equation 

t - tjk [t] = \rj [t) - Tk {tjk {t))\ for aU t < e. 

As follows from Theorem (j3.ip on uniqueness of the retarded times, the function 
tjk is uniquely determined by the above condition. 

Put Tjk — t ~ tjk and call this function the delayed time. Introduce a vector 
function Vjk by the formula 

rjk{t) = Tjit) - rk{tjk{t)) for aU t < e, j ^ k. 

and the unit vector function ejk by 

ejk{t) ^ rjkit)/\rjkit)\ for all t < e, j ^ k. 

Since Tjk — \rj{t) — rk{tjk{t))\ and by assumption the bodies do not collide for 
any t G I the unit vector ejk{t) is well defined. Let Vj{t) be the velocity of the j-th 
body and let yj{t) denote the state vector {rj{t),Vj{t)) and let y{t) be the list of 
all state vectors yj{t) for j = 1, . . . , n. 

Thus the vector function y{t) has its values in the space i?^". By the norm of the 
state vector y G i?^" of the system {rj,Vj), where j 1, . . . , n, we shall understand 

|y| max{|rj|, \vj\ : j = l,...,n}, 

where \rj\ and \vj\ denote the usual Euclidian norm in R^. Since the vector r^fe, 

and also the vector Cjk, is directed from the k-th body towards the j-th body the 
Lorentz force induced by the fc-th body is given by the formula 

qj[Ejk + Vj X {cjk X Ejk)] 

and Newton's force is given by 

—rriQkEjk 

from the fact that force of action and force of reaction are equal in magnitude but 
have opposite directions. Thus the total force Ejk acting onto the j-th body by the 
k-th body is given by the formula 

Fjk = qj[Ejk + Vj X {ejk X Ejk)] - mokEjk 
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Introduce an operator H by the the formula 

(12.1) Hjk{h) = qj{h + Vj x {cjk x h)} — mokh for all h G . 

For fixed y and t and j ^ k the operator Hj^. is linear with respect to h. Denote by 
U the algebra of all linear transformations from the space into itself. 

The total force acting onto the j-th body now can be represented in the form 

(12.2) Fj = ^lll^^. HjkiEjk) + Fj 0, 

where Fjq represents external forces due to some external fields. It could be a 
force field due to self action or due to loss of energy through the electromagnetic 
radiation. We shall want to find conditions that the fields should satisfy in order 
to have uniqueness of the solution of the system of evolution equations. 

Now for a body of rest mass mo moving under the influence of a force F from 
Newton-Einstein formula, time derivative of the momentum is equal to force, we 
must have 

F 

where 7 = (1 — liij^)"^/^ and p = nio^v is the relativistic momentum. Computing 
the derivative of p with respect to time we get 

p = mo'yv + mo'yv = mo'yv + mo7^ {v, v) v 

= mo'y{v + 7^ {v, i)) v) = m{v + 7^ {v, i)) v) 

= mr{v){v), 

where TO ~ mo7 is the relativistic mass and T{v) is the linear transformation of R^ 
into i?"^ given by the formula 

r{v){h) = h + j"^ {v,h)v for all heR^. 

Notice that for fixed velocity v the transformation r{v) represents a symmetric, 
positive define transformation with eigenvalues equal respectively to 

(1 + 7'l^n, 1, 1. 

Thus the inverse transformation r{v) exists and is also symmetric. Since its eigen- 
values are 

(l+7'H^)-\ 1, 1 

and the norm of positive definite symmetric transformation is its maximal eigen- 
value, we must have 

\t{v)\ = 1 for all veR^ 

and we have also 

f(u)r(t;) = e for all v G R^ 

where e denotes here the identity transformation in the algebra Lin{R^,R^) that 
is e{h) = h for all h e R^. 

We shall call the function t; 1-^ F the reciprocal to the function T{v), since 

t{v) = [T{v)]-^ for all v € R^. 
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Thus the Newton-Einstein equation for the j-th body can be written as 

(12.3) r(wj)(wj) = —Fj for ah j = l,...,n 

rrij 

Since trajectories yj = {rj,Vj) are admissible (I2.ip . for stopping time ii = e we 
must have 

sup{|wj(i)| : i < e, j = 1,2, . . . ,n} (? < 1, 

(12 4) 

sup {\vj{t)\ : t < e, j ^ 1,2, ... ,n} =^ A < oo 



Now put x{t) = y{t) for t < s. Such trajectory will be called an initial trajectory. 
Clearly x is of class C^. We shall follow here the notation used in Bogdan [1 . 

Definition 12.1 (Initial domain). For a given initial trajectory x and nonnegative 
numbers q < 1 and A < oo and an interval I = [s,e] define the set 

D^D{x,q,A,I) 

of functions y{t) — {r{t),v{t)), where v(t) — f{t), extending the trajectory x to the 
interval I and satisfying the inequalities 

\rAt)-rJi)\< q\t-t\ for all t, i < e. 
(12.5) ' jyji-'ii I ^ 

\vj{t) ~Vj{t)\ < A\t~t\ for all t, t < e. 

Such a set D will be called an initial domain generated by the parameters x, q. A, I. 

Definition 12.2 (Space C{I,Y) of continuous functions). For any interval I and 
any Banach space Y let C{I, Y) denote the set of all continuous functions f from 
I into Y. 

In the case when the interval / is closed and bounded the space C{I, Y) with 
norm ||/|| = sup{|/(i)| : t & 1} forms a Banach space. 

The set D{x, q. A, I) can be considered as a subset of the Banach space C(I, i?^"). 
It is clear that the original trajectory y of the entire system of bodies belongs to 
the initial domain D{x, q, A, I) for parameters q and A as defined in (|12.4p . 



Definition 12.3 (Space T of admissible trajectories). For fixed n and an interval 
I of the form I — (— cxd, s] denote by T{I) the set of all admissible trajectories x of 
the form 

Xj{t) — (rj(t), Vj{t)) for t £ I and j — 1, . . . ,n. 

Trajectories x of this form will be called admissible trajectories on interval I 
and the set T{I) the space of admissible trajectories on I. 

The union T of all such spaces T{I) will be called the space of admissible 
trajectories, that is 

r-{jT{i). 



Notice that every initial domain D[x, q, A, [s, e]) forms a subset of the space 

r((-cx3,e])cr. 
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Definition 12.4 {SpaceU{Y) of nonanticipating operators) . For fixed n and a fixed 
Banach space Y let P denote an operator on the space T mapping an admissible 
trajectory x on an interval I — {—oo,d] into a continuous function from the space 
C{I,Y). 

Such an operator is well defined if it has the following property: for every two 
functions y, z and every number d d R the condition that y{t) = z{t) for all t < d, 
implies that the images P{y){t) and P{z){t) coincide for all t < d. Every such 
operator will be called a nonanticipating operator and the space U{Y) the space 
of nonanticipating operators with values in the space Y. 

The following theorem represents Th. 4.6 of Bogdan [2]. 

Theorem 12.5 (Initial domain is compact and convex). Every initial domain 
D{x,q, A, I) is compact and convex in the Banach space C — C{I , R^"). 

Thus every initial domain forms a closed set in the space C{I, i?^"). 

Definition 12.6 (Section of Initial Domain). Let D = D{x,q, A, I) be an 
initial domain with I = [s,e] For any d £ [s, e] denote by Ij^ the interval [s, e] and 
by Del the initial domain D{x,q, A, Id). Any such initial domain Dd will be called 
a section of the initial domain D. 

It is plain from the definition that every nonanticipating operator P on an initial 
domain D induces a nonanticipating operator Pd on any section Dd of the initial 
domain D. 

Definition 12.7 (Uniformly Lipschitzian Operator). A nonanticipating operator 
P from an initial domain D = D{x^ q, A, I) into the space C{I,U), of continuous 
functions from the interval I into a Banach space U, will be called uniformly 
Lipschitzian if there exists a constant M such that for very d £ I we have 

\\Pd{v)-Pd{yn\\d<M\\y~y^\\d for all y,y-&Dd, 
where \\ \\^ denotes the norm in the space C{Id,U). 

Obviously the operators y rj and y — >■ Vj are nonanticipating and uniformly 
Lipschitzian on any initial domain D. 

Definition 12.8 (The space Lip{D,U)). Let D = D{x,q, A, I) denote an initial 
domain and U a Banach space. Let Lip{D, U) denote the space of all nonantici- 
pating uniformly Lipschitzian operators from the domain D into the Banach space 
C{I, U). 

If P is an operator in Lip{D, U), let 
\\P\\ =snp{\\P{y)\\ -.yeD}, 

l{P)=mi{M : \\P{y)^P{y~)\\^<M\\y-y-\\^ for all y,y~eD,deI}. 

Clearly since D is compact, this norm \\P\\ is well defined. The norm in the ex- 
pression \\P{y)\\ is understood as the supremum norm in the space C{I, U). Define 

ll^ll +^(^) for all PeLtp{D,U). 
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The space Lip{D,U) with norm \\ ||; will be called the space of nonanticipating 
uniformly Lipschitzian operators. 



13. Nonanticipating differential equations 

Let D = D{x, q, A, I) be an initial domain and A e Lip{D , R^"') . Consider the 
system of differential equations 

(13.1) Vj{t) = Aj{y){t) for all t e I, j ^ 1, . . . ,n. 

Introduce the operator X who's component Xj corresponding to the j-th body is 
defined by the formulas 

Xj{y){t) = {vj{t),Aj{y){t)) for all y ^ {r,v) e D and t e /. 

The right side of the equation p3.1l) represents a continuous function of the variable 
t, so derivative on the left side also represents a contuous function. Thus we can 
apply the integral to both sides to get an equivalent equation 

Vj{t) = Vj{s) + J Aj{y){u)du for all t G /, j = 1, . . . , n. 
We also have 

rj{t) = rj{s) + J Vj{u) du for all t £ I , j — 1, . . . ,n. 
Now if we define the operator ftj for all y E D by the formula 
(rj(s)+ / Vj{u)du, Vj{s) + / Aj{y){u)du) 



we can rewrite the previous two equations in the form 

Vjit) = irj{t),vj{t)) = nj{y)it) for all t e I , j = 1, . . . ,n. 
or even in a shorter form as 

y = ^(y) 

where Q represents the list of operators flj, (j ~ 1, . . . , n). This operator belongs 
to the space Lip{D,R^'^). In this way we reduced the problem to a fixed point 
problem. 

Now we are ready to prove that the fixed point problem has a unique local 
solution. Compare the following theorem with Th. 7 of Bogdan [1]. 

Theorem 13.1 (Local existence and uniqueness). Assume that the operator A is 
nonanticipating and for every value of A > there exists an interval I — [s,e] 
such that A restricted to the initial domain D = D{x,q, A, I) yields an operator 
belonging to Lip{D,R^^^). 

Let Aq > |A(a;)(s)| . Let A > 2Aq. Compactness of the initial domain D = 
D{x,q, A, I) implies the compactness of its image A(D). Thus there exists a non- 
decreasing function uj{6), called modulus of continuity, such that uj(6) — )■ when 
6-^0 and 

\A{y){t)~A{y){t')\<ui\t-t'\) for all y e D, t,t' e L 

Let (3 be such that w(/3) < Aq and X > max{Z(Aj) : j — l,...^n}, where l{Aj) 
denotes the Lipschitz constant of the operator Aj e Lip{D, R^). 
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Let J = [s,s + S] where 6 is such that 

0<6 < min{l,^,A-i} . 

Then the transformation Q maps the domain D = D{x, q, A, J) into D and forms a 
contraction map. Thus for the initial trajectory x of class there exists a unique 
local solution in the domain D to the evolution equation of the system of n bodies 

y{t)=X{y){t) for all tel. 

Proof. First let us establish that the operator O maps the initial domain D into 
itself. Take any y = {r,v) G D and let y = {f,v) = ^{y)- For the position 
component we have 



f{t) = r{s) + / v{u) du for all t < s + S. 

J s 

and for the velocity component 

v{t) = v{s) + / r{u) du for all t < s 

J s 

v{t) = v{s) + / A(y)(u) du for all s <t< s- 

J s 



and 



S. 



Notice that the above formulas represent a trajectory y extending the initial tra- 
jectory X. 

To check the Lipschitz conditions notice that we have for the position component 
\r{t) - f = I / v{u) du\ < q\t - t'^\ for all t, i~ e /. 



For the velocity component we have 

ft 

\v{t) - v{t-)\ 



du 



< 



< 



A{y){s)du 



Hy)is) 



A{y){u) du 



{A{y){s)-A{y){u))du 



A{y){s)du - 

< \A{x){s)\\t - r I + uj{5)\t - r I 

< {Ao + u{d))\t - r I < (^0 + ^{P))\t - r I 

<2Ao\t-t'^\< Alt-f^l for all t,^^ € I, 
Thus we have that y G D. 



The transformation is a contraction. Indeed, take any y^y"" & D and consider 
their images z = fl{y) and z"" = Q{y"'). Introducing the notation for components 
of 2; = {ri,vi) and z~ = {ri,Vi ) we get 

rt rt 



\ri{t)-r^{t)\ = 



/{v{u) — v'"{u)) du = / llw — 
J s 



du 



^ ^\\y — y'" II for all 2/) € D, and t & I 
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[\Aiy)iu)-A{y-){u)) 

J s 



and for velocity component of z, 

\vi{t) - v^{t)\ = l\h(v)(u) - h(v-){u)) du 

< 6X\\y — y'" II for all y, y'" e D, and t G I 
The above inequalities imply 

\\z - z'^W <ma.x{S,6X}\\y - y'^W for all y,y~GD, 
or equivalently 

my)-n{y-)\\<q\\y-y~\\ for all y,y- & D, 

where q = max {S, 5X} . Since it follows from the assumptions of the theorem that 
the constant q is less then 1, the transformation f2 represents a contraction mapping. 

Hence there exists a unique trajectory y G D such that y = or equivalently 
that 

y{t) = X{y){t) foraU tel. 
This completes the proof of the theorem. □ 



14. Existence and Uniqueness of Global Solutions 

Definition 14.1 (Operator with local uniqueness property). Assume that we have 
a nonanticipating operator X -.T ^ i?^" with the property that for every admissible 
initial trajectory x G T there exists a unique local solution of the equation 

y{t)^X{y){t) for all tel 

extending the trajectory x. We shall say that such operator has the local unique- 
ness property. 

Assume that we found a solution y of the differential equation in an initial domain 
D{x,q,A,I), where I = [s,e]. Taking the endpoint e as our new starting point s 
and the trajectory y as our new initial trajectory x one can extend the solution 
onto a larger interval. The object is to show that between all possible extensions 
there exists a maximal extension. Moreover that the maximal extension is unique. 
Here we will make the first step in this direction. 

Theorem 14.2 (Uniqueness of extensions of solutions). Assume that a nonantici- 
pating operator X -.T ^ i?^" has a local uniqueness property for any starting time 

s and any initial trajectory x € 7~((— oo,s]). Let I be an interval starting at the 
point s. On the right the interval may be either open or closed, either bounded or 
unbounded. 

If y and y"" are two solutions of the differential equation 

y{t) = X{y){t) for all t e I, 

extending the initial trajectory x onto the interval I, then these trajectories are 
identical on the entire interval I that is 



y{t)=y-{t) for all tel. 
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Proof. Let / be an interval starting at the point s. On the right the interval may 
be either open or closed, either bounded or unbounded. 

Since difference of two continuous functions on the set / yields a continuous 
function, the function (/) = y{t) — y'^{t) is continuous on /. Consider the set J 
defined by 

J = {i e / : (j){u) = for all u < t.} 

Clearly ,s G J so the set J is nonempty. 

Since / is an interval and since the set J with any two points ti,t2 S J such 
that ti < t2 contains all points u such that ti < u < t2, the set J itself forms an 
interval. 

If the interval J is, either unbounded on the right, or is bounded on the right 
but the least upper bound sup( J) coincides with the right end of the interval I, we 
immediately sec that I = J. 

So let us consider the case when J is bounded on the right but sq = sup( J) does 
not coincide with the right end of the interval I. We shall prove that this leads to 
a contradiction. 

It follows from the continuity of the solutions yi and y^ that they coincide up 
to time t = sq including t = so and are of class C^. Prom our assumption about 
the point sq follows that it must lie inside of the interval /. Taking the part of the 
trajectories up to time t = sq as & new initial trajectory x, we can construct two 
unique local solutions defined on some intervals I\, I2 C / to the right of s. We 
may assume that /i = I2, otherwise we would take their intersection = /i fl I2 
as a new interval. 

Let (5 > denote the length of the interval /'". Since the solution of the equation 
is unique, we would have that 

y{t) = y'^{t) for all t G Iq = I n [so, so + 6), 

which would yield that sq + (5 < so- A contradiction. Thus the interval J must 
coincide with the entire interval /. □ 



Definition 14.3 (Maximal global solution). Assume that X is a nonanticipating 
operator with uniqueness property. If between all possible solutions extending an 
initial trajectory x, there exists a solution with longest possible interval, then such 
a solution is called maximal global solution. 



Theorem 14.4 (Existence and uniqueness of the maximal global solution). Assume 
that X : r is a nonanticipating operator with uniqueness property. 

Then for every initial trajectory x & T, between all possible solutions y of the 
differential equation 

y = x{y), 

which extend x, there exists a maximal global solution. Moreover this solution is 
unique. 
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Proof. To prove the theorem let y denote a trajectory, extending the initial trajec- 
tory X, and defined on some interval I — [s,e]. Assume that y represents a solution 
of the differential equation 

y^x{y). 

Take all such solutions with their intervals / and let / = IJ / be the union of 
all such intervals. If a point t belongs to / it belongs to some interval / being the 
domain of a solution y. Put 

From the theorem on uniqueness of extensions follows that the value y{t) is well 
defined, that is it does not depend on the choice of the function y. 

Indeed if t belongs to any two intervals Ii and I2 then on the intersection J of 
these intervals the corresponding solutions yi and y2 must coincide in accord with 
the Theorem ll4.2l Thus we must have that yi{t) ~ y2{t)- 

Moreover the obtained trajectory y is the solution of our differential equation 
on the entire interval /. Since the graph of the function y contains in it the graph 
of any other solution of our equation, the solution y represents a maximal global 
solution. Again from the Theorem 114.21 follows that this solution is unique. □ 



15. Spaces of nonanticipating uniformly Lipschitzian operators 

We shall remind the reader the definition ll2.8l of the space Lip{D, U) of uniformly 
Lipschitzian operators on the initial domain D. 

Let D = D(x, q, A, I) denote an initial domain and U a Banach space. Let 
Lip{D, U) denote the space of all nonanticipating uniformly Lipschitzian operators 
from the domain D into the Banach space C{I, U), of all continuous functions from 
the interval I into the space U. If P is an operator in Lip{D, U), let 

\\P\\ =sup{||F(y)|| : yeZ?}, 

l{P)^\ni{M:\\P{y)~P{y-)\\^<M\\y-y-\\^ for all y,y- E Da, d E 1} . 

Clearly, since D is compact, this norm \\P\\ is well defined. The norm in the 
expression |l-P(y)|| is understood as the supremum norm in the space C{I, U). Define 

= ll^'ll + ^(^) for all PeLip{D,U). 

The space Lip{D, U) with norm \\ ||; is called the space of nonanticipating uni- 
formly Lipschitzian operators. 

Proposition 15.1 (The pair {Lip{D,U), \\ ||;) forms a Banach space). For every 
initial domain D and for every Banach space U the space Lip{D,U), of nonan- 
ticipating uniformly Lipschitzian operators, equipped with the norm \\ ||j forms a 
Banach space. 

Proof. The proof follows from Th. 2.15 of Bogdan [2]. □ 
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It follows from the above theorem that the space Lip{D, U) is closed under 
addition and scalar multiplication operations induced by respective operations in 
the Banach space U. In the sequel we will extend this class of operations to include 
all Lipschitzian functions. 

The following is Th. 6.10 of Bogdan [2]. 

Theorem 15.2 (Joint Continuity). // P G Lip{D,U) is a nonanticipating uni- 
formly Lipschitzian operator then the map 

{t,y)^P{y){t) 

is jointly continuous from the product I x D into the Banach space U, that is for 
every point (to, yo) € I x D and for every e > there exists S > such that 

\P{y){t) - Piyo)ito)\ < e when \t-to\ < 6 and \\y - yo\\ < S and it,y) e I x D. 



Corollary 15.3 (Uniform continuity of P{y){t)). Assume that P G Lip{D,U). 
Since the set I x D, as product of compact sets, is itself compact. Thus the function 
{t, y) — >■ P(y)(t) is uniformly continuous on the set I x D into the Banach space U. 

Definition 15.4 (The set ran{P) is the range of operator P). Assume that P e 
Lip{D,U). For shorthand we shall write ran{P) to denote the set 

P{D){I)^{P{y){t):yeD,teI}. 

Corollary 15.5 (Compactness of ran{P)). Assume that P e Lip{D, U). Since the 
set I X D is compact and the function {t,y) — > P(y)(t) is continuous on it, the 
set ran{P), as an image of a compact set by means of a continuous map, forms a 
compact set in the Banach space U. 



Definition 15.6 (Lipschitzian Operators over Cartesian Products). Let 

Ui, . ■ . ,Un\ Z 

be some Banach spaces. We shall say that an operator P, from a subset G of the 
cartesian product Ui x . . . x Un into the Banach space Z, is Lipschitzian on G, if 
for some constant M we have 

\\P{UI, ...,Un)- P{U^, . . . ,0|| < M(||U1 - wrll + • • • + \\Un " <||) 

for all (ui, . . .,Un), (u^, ■ ■ ■,«;;') S G. 



The following is Th. 6.15 of Bogdan [2]. 
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Theorem 15.7 (Lip(D, U) spaces arc closed under composition with Lipschitzian 
operators). Let D be an initial domain and Ui, . . . ,Un; Z be some Banach spaces. 
Assume that Pj G Lip{D,Uj) for j = l,...,n and Q is a Lipschitzian operator 
from the product ran{Pi) x • • • x ran{Pn) into the space Z. Define the composed 
operator P by the formula 

P{y){t) = Q{Pi{ym, Pn{y){t)) for all yeD and t e /. 

Then the operator P belongs to the space Lip{D,Z). We shall use a shorthand 
notation for the operator P = Q o (P^, . . . , P„) or shorter P = Q{Pi, ■ ■ ■ , Pn) when 
this does not lead to a confusion. 

Definition 15.8 (n-linear bounded operators ). Let Ui, . . . , Un, Z be some Banach 

spaces. We shall say that an operator P, from the Cartesian product Ui x . . . x Un 
into the Banach space Z, is n-linear if for every variable uj, where j = 1,. . . ,n, 
when other variables are fixed the map 

Uj P{ui,...,Uj,...,Un) 

is linear from the space Uj into the space Z. Such an operator is said to be bounded 

if for some constant M we have 

|P(wi, . . . < M |u2|, • • • , l^nl for all Ui £ Ui, . . . ,Un & Un. 

Proposition 15.9 (n-linear bounded operator is Lipschitzian on bounded sets). 
If P is an n-linear bounded operator from a Cartesian product Ui x ■■ ■ x Un of 
Banach spaces into a Banach space Z, then it is Lipschitzian on every set of the 
form 

Bs = {{u-L, ...,Un)eUiX---xUn: \uj\ < 5, j = l,...,n} 
Proof. The proof is straightforward and we leave it to the reader. □ 

Corollary 15.10 (Compositions with n-linear operators). // Q is an n-linear 

bounded operator from a Cartesian product Ui x ■ ■ ■ x Un of Banach spaces into 
a Banach space Z, and Pj are in Lip(D, Uj)for j = 1, . . . ,n, then the composition 
Q o (Pi, . . . , P„) belongs to the space Lip{D, Z). 

Proof. Since the real- valued functions Uj \uj\ arc Lipschitzian on Uj, they are 
continuous. Thus on the compact set ran{Pj) they attain their maximum 5j. Put 
5 = max {6j : j — \, . . . ,n\ . 

Then the operator Q is Lipschitzian on the bounded set Bs and the operator 
representing the composition of operators 

P = Qo(Pi,...,P„) 

is well defined and we have P G Lip{D, Z). 

□ 
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Definition 15.11 (DifFerentiable operators). Let f be an operator from an open 
set G in a Banach space U into a Banach space Z. We shall say that the operator 
f is dijferentiable at a point a £ G if there is a ball 

B{a, r) = {u eU : \u- a\ <r} <Z G 

and a linear bounded operator g e Lin{U, Z) such that 

\f{x)-f{a)-g{x-a)\<o{\x-a\) for all x € B{a,r), 

where o{h) denotes some function of the variable h > such that 

o{h)/h-^0 when — > 0. 

Such an operator g is unique, and is called the derivative of f at the point a, and 
it is denoted by f'{a) = g. 

It follows from the above definition that if the map / has a derivative at a point 
a € G then / is continuous at a. Notice also that the space Lin{R, R) is isomorphic 
and isometric with the space R of reals. 

The theorem, that follows, can be found in Cartan "W, page 27, Th. 2.2.1]. 

Theorem 15.12 (Chain Rule). Assume that U, V, W are Banach spaces and G C 
U and H <Z V are open sets. Let f : G ^ H and g : H ^ W be differentiable on 
their domains. The composed function h = f o g defined by 

h{u) = / ° 9{u) = f{g{u)) for all ueG 

is differentiable and 

h'{u)^f'{g{u))g'{u) for all u e G. 
The following theorem can be found in Cartan [HJ page 41, Th. 3.3.2]. 

Theorem 15.13 (Maps with bounded derivative on convex sets are Lipschitzian). 

Let f be an operator from an open set G in a Banach space U into a Banach space 
Z. Assume that the derivative f'{u) exists at every point u ^ G. Assume that for 
some convex set W <Z G and a constant M we have 

\f'{u)\<M for all uC^W. 

Then 

\f{u)~ f{u~)\<M\u~u'-^\ for all u,u~eVK. 

The space Lin{U, U) of linear bounded operators from a Banach space U into 
itself beside being closed under addition and scalar multiplication is also closed 
under the composition of operators: P o Q e Lin{U, U) for all P,Q^ Lin{U, U). 
This operation has the property that l^'oQI^I-PllQI. 
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Definition 15.14 (Banach Algebras). A Banach space U is called a Banach al- 
gebra, if it is equipped with a bilinear operation {u,w) uw, from the product 
U X U into U, that is associative: (uw) z = u (w z), and such that 

\uw\ < \u\ \w\ for all u, w €z U. 

If in addition there is an element e in U such that en = ue = u for all u d U, 
then such an algebra is called a Banach algebra with unit. In such a case an 
element u is called invertible if for some w €z U, called the inverse of u, we have 

uw ^ w u = e. 

The unit element and the inverse elements are unique. We denote the inverse of u 
by u~^. 

The following proposition is very useful in numerical computations of the inverse 
transformations. It represents a simple application of the Banach fixed point the- 
orem and it is essential in the development of the theory presented in this paper. 

It represents Pro. 6.25 of Bogdan [2]. 

Proposition 15.15 (Inverse {e — v)~^ exist if |u| < 1). Let U be a Banach algebra 
with unit. Then for every element v € U such that \v\ < 1 the inverse w = {e — v)~^ 
exists. It is a fixed point of the operator f defined by the following formula 

f{u) ^ e + vu for all u (z U. 

Moreover the sequence Wn — where Wq = is of the form 

w„ e + V + v'^ -\ + = v'' for all n> 

and we have the following estimate for the distance of the fixed point w and the 
approximation Wn 

\w — Wn\ < TT 0,11 n > 0, 

1 - \v\ 

and we also have an explicit formula for the inverse element w as the sum of an 
absolutely convergent series 

w = e + V + + ■ ■ ■ — w". 

Notice that in any Banach algebra Lin{U, U) the identity map; e(u) = u for all 
u € U, is the unit element. It is easy to prove that the norm of the unit element e 
in any Banach algebra is equal to 1, |e| = 1. Therefore for every invertible element 
u we have 1 = < |m| thus \u^^\ > 0. 

The following is Pro. 6.26 of Bogdan [2]. 

Proposition 15.16 (When the inverse map u i-> is Lipschitzian?). Assume 
that U is a Banach algebra with unit. 

(1) If for some uq the inverse u^^ exists then every element in the ball 
B{uq, r) ^ {u eU : |u - wo| < r} , 
where r — 1/\uq^\, has an inverse. 
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(2) The dom,ain G of existence of the inverse is an open set and the function 
f(u) ~ u^^ is continuous on G. 

(3) The function f{u) = is differentiable on G and 

f'{u){h) = —u~^hu~^ for all h £ U and u £ G. 

(4) If W is a compact convex set such that W C G, then f{u) is Lipschitzian 
on W. 

Theorem 15.17 (If reciprocal g oi f G Lip{D,U) exists then g £ Lip{D,U)). 
Let D be any initial domain and U a Banach algebra with unit e. Then if f is a 
uniformly Lipschitzian function from the domain D into the algebra U and 

\3{y){t)] [f{y){t)] = e for all yeD andtel 
then g G Lip{D, U). 

Proof. By continuity of the functions (y^t) i— > f{y){t), and x i— > x~^, and the norm 
y \y\ and compactness of the set D x I we get the existence of a constant M 
such that 

I [fiy){t)]~^ \<M for all ye D and t e /. 
From the algebraic identity 

we get 

\9{yim - giy^mi = |[/(2/i)(<)]-M/(2/2)W - /(2/i)W][/(2/2)(t)]-^| 

< M^l{f) \\yi - yalid for all yi, y2 e and t € d € I 

Therefore 

l{g) < MH{f) 

and this implies that g e Lip{D, U). □ 

Theorem 15.18 (/ G Lip{D,B!') ^ fj G Lip{D,R) V j = l,...,k). Let f be 
a map from an initial domain D into the space R''. Assume that / = (/i, . . . , fk) 
where fj is from D into the reals R for all j — 1, ... ^k. 

Then f belongs to the space Lip{D, R^) if and only if each component fj belongs 
to the space Lip{D,R). 

Proof. The proof is obvious and we leave it to the reader. □ 

Theorem 15.19 (/ G C\R'',R"') and ,g G Lip{D,R^) => f o g e Lip{D,R"')). 
Assume that f is a continuous function from the space R^ into the space R™ having 
continuous partial derivatives of order 1. 

Then if g is a uniformly Lipschitzian function from an initial domain D into the 
space R^ then the composed function f o g is uniformly Lipschitzian from D into 
the space R^. 
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Proof. By assumption the derivative f'{x) exists and is continuous from R'' into 
the space Lin{R^,R™') of Hnear transformations from R^ into the space i?™. Thus 
on every closed bounded convex set the map x >■ is bounded. Thus from 

Th. 3.3.2, page 41, of Cartan [6 follows that / is Lipschitzian on every such set. 

Therefore from Th. 1 15. 71 follows that f og is uniformly Lipschitzian on the initial 
domain D. □ 



As a consequence of the above theorem we can conclude that for any polynomial 
p the following is true. 

Corollary 15.20 {xj G Lip{D, i?) V j => p{xi, . . . , Xk) E Lip{D, R)). Assume that 
D is an initial domain. Assume that p{x) is a polynomial of coordinates Xj of the 
vector X. If components Xj = Pj G Lip{D, R) then p o (xi, . . . , Xk) is in the space 
Lip{D,R). 

16. Initial Domains with Positive Separation 

In this section we will derive some estimates involving the speed c of light, so 
for a while we will explicitly use this constant. Only in later sections it will be 
more convenient to change the units. Compare the following developments with 
constructions in Bogdan 1 in particular with results in Th. 10 and Th. 11. 

Let D = D{x, q, A, I) be an initial domain. In this section we will develop tools 
to study further the properties of the spaces Lip{D, U) that are need for analysis of 
operators generated by moving bodies in gravitational and electromagnetic fields. 

Definition 16.1 (Lattice Operations). Assume that R as before denotes the field 
of real numbers. Introduce operations V and A from R x R into R by the formulas 

p\j r ~ max {p, r'\ ~ -^{p + r + \p — r\), p Ar ^ niin {p, r} = —{p + r— \p — r\). 

Any linear space U of real-valued functions closed under composition with these two 
operations will be called a linear lattice. 

Notice that if f £ U then the absolute value \ f\ of the the function f is also in 
U, since \ f\ = (— /) V /. Now, if in addition the space U is a Banach space with a 
norm \\ || , satisfying the following implication 

5= I/I =^ llffll < for all feu, 

then U is called a Banach linear lattice. It is easy to see that the last condition 
guarantees that the lattice operations are continuous in the norm topology. 

The operations V and A are commutative, associative, and Lipschitzian, since 
the unary operation p — > |p| is Lipschitzian. The following theorem is useful in 
getting estimates of Lipschitz constant for operators in various spaces Lip{D, U). 

The following represents Th. 7.2 of Bogdan [2]. 
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Theorem 16.2 {Lip{D,R) is a Banach linear lattice and a Banach algebra). For 

every initial domain D the Banach space Lip{D,R), with the norm defined by 

||P||;= ||P|| + for all PeLip{D,R), 

where 

||P||=sup{||P(2;)|| -.yeD}, 

l{P)^mi{M : \\P{y)~P{y-)\\^<M\\y-y-\\^ for all y,y~eDd,deI}, 
forms a Banach linear lattice and a Banach algebra with unit. 

In the sequel we will need a separation and minimal time delay operators. 
Definition 16.3 (Separation and Minimal Time Delay Operators). Assume that 

D = D{x,q,A,I) 

represents an initial domain with I = [s,e]. Define the operator sep, called the 
separation operator, by the formula 

sep(y)(i) = min{|rj(t) - rfe(t)| : fc = 1, . . . ,n, j 7^ fc} for all yeD,t<e, 
and the operator mtd, called the minimal time delay operator, by the formula 
mtd{y){t) = mm{Tjk{t) : j,k ~ 1, . . . ,n,j ^ k} for all y ^ D, t < e. 

The following represents Pro. 7.4 of Bogdan [2]. 

Proposition 16.4 (Operators sep, mtd belong to Lip{D,R)). The separation 
operators sep and mtd, when the trajectories y E D are restricted to interval I, 
belong to the space Lip{D, R) of nonanticipating uniformly Lipschitzian operators. 

Proof. The proof follows from theorem 116.21 □ 

Definition 16.5 (Time Delay Bound and Separation Bound). Let D — D{x, q, A, I) 
be an initial domain. The time delay bound th for the domain D is defined by 

tb = inf {mtd(y)(i) : y e D, t e 1} . 

Similarly we define the separation bound Sb for the domain D by 

Sb = inf {aep{y){t) : y£D,teI}. 

Notice since ranges of the operators mtd and sep are compact the constants tb 
and Sb are well defined and are finite. 

Notice that tb = if and only if Sf, — 0. Indeed since the map (y, t) — mtd{y){t) 
is continuous and the set D x I is compact, there exist a trajectory y £ D and a 
point t € I and a pair of indexes j, k such that for some j, k we have 

O^tb^ Tjk{t) = -^\rj{t)-rk{t-Tjk{t))\ = \j{t)-rk{t)\ > hep{y){t) > ^Sb > 0, 
thus Sb — 0. The converse that Sb — implies tb = can be proved similarly. 

Definition 16.6 (Domains with Positive Separation). Any initial domain D such 
that Sb > will be called a domain with positive separation. 

The following represents Pro. 7.7 of Bogdan [5]. 
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Proposition 16.7. If D ^ D{x,q, A, I) is an initial domain with positive separa- 
tion and I = [s,e], then for every d such that < d — s < tb and d < e the emission 
time operators tjk satisfy the inequality 

tjk{t) ~ t — Tjk{t) < 8 for all y G and t G [s, d]. 

Proof. Notice that for the section Dd of the domain D the time interval is [s,d]. 
Thus for any t G [s, d] and any trajectory y G we have 

tjk{t) = t - Tjk{t) < d-tb < d - {d - s) ^ s. 

This completes the proof. □ 

The operators y rk{tjk), and y -J> Vk{tjk), and y ~^ Vk{tjk), are constant on 
the domain Dd and they are uniquely determined by the initial trajectory x, since 
tjkit) < s and by definition the initial trajectory x has a continuous derivative. 
Therefore each, of the above operators, represents a continuous function of the 
variable t. Thus all the above operators belong to the space Lip{Dd, R^)- 

Notice that the time delay bound t^ for the domain D and the time delay bound 
t'lj for its section Dd satisfy the inequality tb < t'^^. Thus taking eventually the section 
Dd as our new domain we observe that the following corollary must be true. 

Corollary 16.8 (Operators rfe(ijfe), Vk{tjk), v'^.{tjk), r^k are in Lzp(D, i?^)). As- 
sume that D = D{x,q, A, I) is an initial domain with positive separation, and with 
I — [s, e]. If we have that < e — s < tf,, then each of the following operators 

rk{tjk), Vk{tjk), v'kitjk), rjk ^ (rj -rkitjk)), (vj ~ Vk{tjk)), 
where j,k — 1, . . . ,n, j ^ k, belongs to the space Lip{D, R^). 

The following represents Th. 7.9 of Bogdan [51. 

Theorem 16.9 (sb < c(l + q)tb)- Assume that D = D{x, q, A, I) is an initial 
domain, and Sb is the separation bound for the domain D, and tb is the time delay 
bound for the domain D. Then 

Sb < c(l + q)tb. 

The following represents Th. 7.10 of Bogdan [2]. 

Theorem 16.10 (Lower estimate of separation bound s^). Let D ~ D{x, q, A, I) be 
an initial domain with I = [s,e] and with positive initial separation sep(a;)(s) > 0. 
// the length of the interval I ^ e ~ s is such that 

sep(x)(s) 

e- s < 5 , 

3cq 

then we have the following lower estimate for the separation bound Sb 

sep(a;)(s) 



The following represents Th. 7.11 of Bogdan [2]. 
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Theorem 16.11 (Negative powers of time delays are in Lip{D,R)). Let D = 
D{x, q, A, I) be an initial domain with positive separation. Then for every pair of 
indexes j ^ k the operator y — > T^^", where rt > 0, is in Lip{D,R). Its Lipschitz 
constant can be estimated by 

(16.1) l{T-,n < ntr'-H(T,,) < ^(Y^^- 

The following represents Cor. 7.12 of Bogdan [5]. 

Corollary 16.12 (Operators \rjk\, I^jTcI £jk)- Assume that the initial domain 
D = D[x, q, A, I) is with positive separation. If the length of the interval I — [s, e] 
is such that e — s < ti,, then the operators 

V \rjkl y kj/cT", {n > 0) for all j 

are in Lip{D,R), and the unit operator 

y ^ Cjk = \rjk\~^rjk for all j^k 

is in Lip{D, R^). 

17. LIP(D,U) SPACES OVER REGULAR INITIAL DOMAINS 

Definition 17.1 (Regular initial domain). An initial domain D = D{x,q, A, I), 
where I = [s,e], with positive initial separation sep(x)(s) > and such that 

0<e-.<!2PM(f) 
3(l + g)c 

will be called a reguleir initial domain. 

Since on every regular initial domain D we have the inequalities for the time 
delay bound 

<t,= inf {tdb(y)(i) : yeD,teI} 

3(1 + q)c 
and for the separation bound 

'''^^3^^'^ <s, = inf {sep(y)(i) : yeD^tel}, 

the bounds t^ and si, are positive. It will be convenient from now on to assume, if 
not specified otherwise, that the domain D is regular. 

From now on we will assume that c — 1. We have already proved various parts 
of the following theorem but to be sure that we did not overlooked any particular 
fundamental field we collet them all into one theorem and go over the proofs again 
now that we have more tools to do that. 

Theorem 17.2 (Fundamental fields due to fc-th body acting onto j-th body are 
in Lip[D,R"^) for m = 1 or m = 3). Assume that D = D(x, q, A, I) is a regu- 
lar initial domain. Consider the fields acting onto the j-th body by the k-th body 
corresponding to the fundamental fields of the system of the trajectories. 
Then the scalar fields 

y ^ Tjk, Ujk, tjk, Zjk 
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all belong to the space Lip{D, R) and the vector fields 

y^rj, Vj, rjk=rj~rk{tjk), e^-fe = Uj-fer^fe, rj{tjk), Wj(ijfe), aj(tjfe) 
all belong to the space Lip{D,R^). 

Proof. It is obvious that operators y h- ^ rj and y i— >• Vj are in Lip{D,R^). The 
remote operators, as fonows from ()16.8p . 

V ^ rj{tjk), Vj{tjk), aj{tjk) 
are also in Lip{D, R^). From Hnearity of Lip{D, R^) fohows that operator 

y rjk = rj - rkitjk) 

is in Lip{D,R^). The operator y ^ T = \rj— rk{tj}~)\ as a composition of Lips- 
chitzian function r i— > \r\ with rjk is in Lip{D, R). 
Since on a regular domain we have the estimate 

Tjk{y){t) >tb>0 for all y e D and t e I 

the operator Ujk as the reciprocal operator to Tjk is in Lip{D, R). Thus the operator 
Cjk = Ujkrjk considered as composition of operators Ujk G Lip{D, R) and operator 
rjk S Lip{D,R^) with bilinear bounded operator (A,r) i— > Ar is in Lip{D,R'^). 

Similarly we can prove that operators tjk = t — Tjk and Zjk = (1 — {^jk, i^jfc))"^ 
are in Lip{D, R). □ 



Theorem 17.3 (Newton-Feynman fields Ejk are in Lip{D,R^) ). Assume that 
D = D{x, q, A, I) is a regular initial domain. Consider the Newton-Feynman field 
Ejk acting onto the j-th body by the k-th body of the system y (z D of the trajectories. 
Then for every pair of indexes j ^ k the field Ejk G Lip[D, R^). 



Proof. According to Bogdan-Feynman Theorem the Newton-Feynman field 

E generated by a single trajectory and expressed in terms of the fundamental fields 
is given by the formula 

E = —uz^a + uz^{e, a)e — uz^{e, a)v 

(-^'^■-^) 23 2 3/ \ 23 2 3/ \ 

+ u z e — u z {v,v)e — u z v + u z {v,v)v. 

Thus the components of = {Ei, E2, E3) are representable as polynomials of the 
components of the fundamental fields. Hence Ejk is representable as a polynomial 
of the coordinates of all the fields appearing in the previous theorem. So each 
coordinate of the field Ejk is in Lip{D, R) so the field itself Ejk is in Lip{D, R^). □ 
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18. Existence of solutions to the evolution equations in force fields 
of special theory of relativity 

We remind the reader that the Newton-Einstein (|12.3|) equation for the j-th 
body can be written as 

r{vj){vj) = —Fj for ah j = 1, . . . ,n 

or equivalently as 

(18.1) Vj = —t{vj){Fj) for ah j = l,...,n 

rrij 

Definition 18.1 (Regular nonanticipating operator). By a regular nonantic- 

ipating operator we shall understand an operator F : T ^ R^" such that for 
every regular initial domain D ~ D(x,q, A, I) the restriction of F to D represents 
a uniformly Lipschitzian operator, that is, it represents an element of the space 
Lip{D,R^'^). 

Notice that in the following theorem the nature of the force operator F is not 
important. Important is just the mathematical assumption that it represents a 
regular nonanticipating operator. This is the property we were looking for. 

Tiieorem 18.2 (Maximal solutions of equations of special theory of relativity). 

Assume that the operator F -.T ^ R^" representing the force field forms a regular 
nonanticipating operator. 

Then for every admissible initial trajectory x of the system of n bodies such that 
at the initial time t — s the separation operator satisfies the condition 

sep(a::)(s) > 0, 

there exists a unique maximal solution y to the Newton- Einstein system of equa- 
tions ^ 

Vj ^ t{vj){Fj) for all j^l,...,n 

nij 

satisfying the condition 

sep{y){t) > for all t e I 
and extending the initial trajectory x onto the time interval I . 

Proof. Let U denote the algebra of all linear transformations from into R^ . The 
map 

g{r, u,x) — r u{x) for all r^R, u^U, x^R^ 

is trilinear and bounded. 

Take any regular initial domain D ~ D{x, q, A, I). Consider first the map 

y ^ "17^ = ™oj-^\/o~~kF)- 

Notice that Pj{y) = (1 — \vj\'^) represents a polynomial of the coordinates of the 
vector operator Vj E Lip{D,R^). Thus pj G Lip{D,R). We also have the estimate 
a — (1 — g^) < (1— l^jp) < 1. Since the square root function is differentiable on 
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the closed interval [a, 1] and has continuous derivative, it is Lipschitzian. Thus the 
operator y i— > nij^ is in Lip{D, R) for every D. 
Consider the operator y — > T{vj). By definition 

r{vj ){h)^h + 7^ {vj , h) Vj for all he R^. 

Notice 7^ = pj^ . Since we established before that p e Lip{D,R) its reciprocal 7^ 
is in Lip(D, R). 

Now consider the 4-linear bounded map 

S{r,vi,V2, h) — r {vi, h) V2 for all r £ R; vi,V2,h£R^ 

and define map 

S{r,vi,V2){h) = r {vi, h) V2 for all r e R; vi,V2,heR^ 

Clearly the map S is trilinear and bounded and maps the triple {r,vi,V2) into an 
element of the algebra U. Thus wc have S{'y'^,Vj,Vj) S Lip{D,U). Since the unit 
element e of the algebra [/ as a constant operator belongs to Lip{D, U) and 

T(v.j) = e + S{-/],Vj,Vj) 

we must have T{vj) G Lip{D, U). Therefore the reciprocal operator y i— > T{vj) must 
belong to Lip{D, U) for every regular initial domain D. 

Hence the operator ^r(wj)(_Fj) belongs to the space Lip{D,R^) for every reg- 
ular initial domain D. Therefore the operator y 1— -^i'ivjjiFj) forms a regular 
nonanticipating operator. 

The above on the basis of the theorem on global maximal solution leads to the 
conclusion stated in the theorem: There is a unique maximal solution satisfying 
the separation condition and extending the initial trajectory. □ 



19. Existence of solutions to equations of gravitational 
electrodynamics with external force fields 



We have introduced in ()12.ip operators Hjk by the the formula 
(19.1) Hjk(h) = qj{h + Vj x (e^fe x h)} - niokh for aU h G R^ . 

These operators have the following property. 

Proposition 19.1 (Operators Hjk £ Lip{D,U) V j 7^ k). The operators Hjk 
for every pair of indexes j ^ k are well defined on any regular initial domain 
D = D{x, q, A, I) and belong to the space Lip{D,U) where U denotes the alge- 
bra of all linear operators from R^ into R^. Thus the operators Hjk are regular 
nonanticipating operators. 

Proof. The proof is similar to the proof concerning the operator T{vj) and we leave 
it to the reader. □ 
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Definition 19.2 (Lorcntz-Newton force operator). By the Lorentz- Newton force 
operator we shall understand the system of the operators, for j^l,. ■ ■ ,n, given by 
the formulas 

(19.2) Fj ^Y!''" H,k{Ejk) for all j = l,...,n, 
where 

(19.3) Hjk{h) = qj{h + Vj X (ejk x h)} — mofe/i for all h <= 

and Ejk represents the Newton- Feynman field acting onto the j-th body from the 
trajectory of the k-th body of the system. 

Theorem 19.3 (Equations of gravitational electrodynamics with external force 
field). Assume that Fjq : T >■ R'^" are regular nonanticipating operators represent- 
ing external force fields. Assume that Fj (j — 1, . . . ,n) denotes the Lorentz- Newton 
force operator. 

Then for every admissible initial trajectory x, defined on the time interval (— oo, s], 
with positive separation sep(a;)(s) at time t = s, there exist an interval I staring at 
s, such that the Newton- Einstein system of differential equations 

Pj{t) — Fjit) + Fj o(i) for all t €z I and j — 1, . . . ,n 

has a unique maximal solution having positive separation 

sep{y){t) > for all tel. 

Proof. Let U denote the algebra of all linear operators from R'^ into R^ . Notice that 
Lorentz-Newton force operator forms a regular nonanticipating operator. Indeed, 
for fixed j ^ fc the operator Hji. G Lip{D, U), where D = D{x, q, A, I). 

The operator Ejk G Lip{D,R^). Since the map S{u,r) = u{r) considered on 
the product U x R^ into R^ is bilinear and bounded we have that Hjk{Ejk) G 
Lip{D,R^). 

Now from linearity of the space Lip{D, R'^) follows that the Lorentz-Newton 
force Fj G Lip{D,R^). Since the force FjQ represents a regular nonanticipating 
operator the sum Kj = Fj + _Fjo belongs to Lip{D,R^). Thus the equation (|19.3p 
is equivalent to equation 

^j(^) = [—^jW) foi' all tG I and j = 1, . . . , n. 

ruj 

We have proved in the preceding section that such system of equations has a unique 
maximal solution satisfying the condition that separation between the bodies is 
positive at any time t G I. □ 

20. Computer programs for DERIVE software system 

The results presented in this paper were obtained with the help of an interactive, 
symbolic, DOS based, computer program DERIVE, version 2.06, developed by Soft 
Warehouse Inc. in Honolulu, Hawaii. 

Since the differential formulas obtained in this paper essentially yield polyno- 
mials of the involved field variables, any program that can handle operations on 
polynomials of an arbitrary number of variables can be used. 

For a reader who would be interested in verifying the formulas we enclose the 
header and the actual programs used in the computations. 
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To understand these programs we need to introduce rudimentary elements of 
grammar used by DERIVE. 

<variable> consists of a string of lower case letters a-z, digits 0-9, and underscore 
character _ . A variable must start with a letter. 

<algebraic expression> consists of variables joint by usual operations: addi- 
tion + , subtraction - , multiplication * , division / , and exponentiation ~ . Paren- 
thesis ( ) can be used to indicate the order of operations. 

<assignment statement> has the form 

<variable> : =<algebraic expression><cr> 

where <cr> denotes the carriage return character. 

<comment> consists of any text enclosed into double quotes "<text>". 

Using tilde ~ as the continuation character one can break long assignment state- 
ments or comments between lines. On output DERIVE breaks line after every 80 
characters. 

DERIVE accepts input files containing programs. For illustration assume that 
we prepared using a plain ASCII text editor a file hh.mth containing the text as 
in the program below representing the header with basic formulas. Notice that 
the sequence in which the formulas are entered in the program is important. One 
should not use things that where not defined in previous formulas. 

To make a test to see that DERIVE understands the program run DERIVE. 
On entering derive a window appears with menu commands at the bottom of the 
window. Notice that each command has one letter in upper case. To select such a 
command just press the corresponding letter on the keyboard. 

To load the program select from consecutive menus: Transfer, Load, Derive, 
hh. Press <cr> that is <Enter> on most keyboards. DERIVE will load the file 
hh.mth and the last statement in the file will be highlighted. To browse through 
the statements use directional keys <up> and <dowii> on the keyboard. The selected 
statement will be highlighted. To enter a statement to examine its components use 
directional keys <right> and <lef t> and to enter or exit a subcomponent use the 
keys <up> and <down>. 

On a highlighted component one can perform any operation from the menu. We 
need just one operation: Expand. It permits one to expand a compound component 
into a polynomial of variables that are involved in the component. 

Expanded components are appended at the bottom of the file hh.mth To save 
the computations select from the consecutive menus commands: Transfer, Save, 
Derive, <cr>. The file hh.mth will be saved to its original location overwriting the 
previous version. To exit DERIVE select from the menu command Quit. 

Using your favorite ASCII editor update the file by new statements and com- 
ments. One can repeat this process until the goal of the computations is met. 

In the following programs the names of the variables were preserved, if possible, 
as they appear in the paper with proper adjustments to lower case letters. 

To represent the dot product {e,v) of vectors e and v, we use variable edv To 
represent a component Vi of the vector v we use v_i. 

The prefix d_ in front of a variable represents the partial derivative D and the 
prefix d_i the partial derivative Thus De and Die are translated as d_e and 
d_ie, respectively. 

Header used in the proofs of the results of the paper 
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This header was used in front of all the programs that follow. 
"Header for electromagnetism" 

M M 

"Define basic identities for dot product:" 
ede:=l 

"To avoid double entrees for dot product veiriables, use" 
"notation that follows order of vector variables as in the list 
"e, V, a, dot_a; thus use vda and not adv for dot product of" 
"vector a with vector v. So edv is ok but not vde." 
edv:=l-l/z 
II II 

"Derivatives with respect to D_i" 

d_it : =z*e_i 

d_itau:=-z*e_i 

d_iu:=-u"2*z*e_i 
II II 

"Meaning of variables in the following expressions" 

" un_ii=l; ee_ii=e_i*e_i" 

" ev_ii=v_i*e_i=e_i*v_i" 

d_ie_i : =+u*un_ii-u*z*ee_ii+u*z*ev_ii 

" ea_ii=e_i*a_i" 

d_iv_i : =-z*ea_ii 
II II 

"Notice that edot_a_ii=e_i*a'_i" 

d_ia_i : =-z*edot_a_ii 
II II 

"un_i denotes here the i-th unit vector of standeird base" 
d_ie : =+u*un_i-u*z*e_i*e+u*z*e_i*v 
d_iv:=-z*e_i*a 
d_ia: =-z*e_i*dot_a 

d_iz : =+u*z~2*v_i-u*z*3*e_i+u*z'2*e_i+u*z"3*vdv*e_i-z"3*eda*e_i 

"D_iz:=+uz*2v_i -uz*3e_i +uz"2e_i +uz"3(v,v)e_i -z"3(e ,a)e_i" 
II II 

"Derivatives of dot products:" 

d_iede : =0 

d_iedv : =-u*z*edv*e_i+u*v_i+u*z*vdv*e_i-z*eda*e_i 

d_ieda: =-u*z*eda*e_i+u*a_i+u*z*vda*e_i-z*eddot_a*e_i 

d_ivdv : =-2*z*vda*e_i 

d_ivda:=-z*ada*e_i-z*vddot_a*e_i 

d_iada:=-2*z*addot_a*e_i 
II II 

"Derivatives with respect to time" 

d_tau:=z 

d_t:=l-z 

d_u:=u'2*z-u'2 

d_z:=u*z-2*u*z'2+z'3*eda+u*z'3-u*z~3*vdv 
II II 

"Time derivatives of vector fields" 
d_v:=z*a 

d_e : =-u*e+u*z*e-u*z*v 
d_a:=z*dot_a 
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"Time derivatives for components of the vector fields" 
d_v_i : =z*a_i 

d_e_i : =-u*e_i+u*z*e_i-u*z*v_i 

d_a_i : =z*dot_a_i 
II II 

"Time derivatives of dot products:" 
"Remember the ordering list: e, v, a, dot_a" 
d_ede : =0 

d_edv : =-u*edv+u*z*edv+u*z*vdv+z*eda 

d_eda : =-u*eda+u*z*eda-u*z*vda+z*eddot_a 

d_vdv : =2*z*vda 

d_vda : =z*ada+z*vddot_a 

d_ada : =2*z*addot_a 

"End of header" 

II M 

JJJJJJJ>J>J>J>J>J>J>J>J>J>J>J>>J>J»J»J»»J»JJ»J»»J»J»»J»JJ»J»JJ»J»»J»JJ»J»JJ 



Check that scalar potential (p = uz satisfies the wave equation 

Include the header at the beginning of this file to process the program. 

"Checking that scalar potential Phi=u*z satisfies the wave equation" 
II II 

"Compute time derivative D(u*z)" 
d_phi : =z*d_u+u*d_z 

"Expanding right side of above yields" 

d_phi:=-u"2*z"3*vdv+u"2*z"3-u"2*z*2+u*z"3*eda 
II II 

"Compute derivative D of each term of the above sum" 

ttl:=-(2*u*d_u*z"3*vdv+u"2*3*z~2*d_z*vdv+u"2*z'3*d_vdv) 

tt2 : =2*u*d_u*z*3+u*2*3*z~2*d_z 

tt3 : =- (2*u*d_u*z~2+u*2*2*z*d_z) 

tt4 : =d_u*z*3*eda+u*3*z"2*d_z*eda+u*z*3*d_eda 
II II 

"Thus the second derivative D~2 of u*z is" 

dd.phi : =ttl+tt2+tt3+tt4 
II II 

"'/.'/,'/.'/,'/.'/,'/.'/,'/. Now compute partials D_i of u*v '/X'/X'/X'/XV 

d_iphi : =z*d_iu+u*d_iz 

"Expanding right side of above yields" 

d_iphi :=+u*2*z*3*vdv*e_i-u*2*z*3*e_i+u*2*z"2*v_i-u*z"3*eda*e_i 

II n 

"Partial derivative D_i of each term of the above sum is" 

uul : =+2*u*d_iu*z*3*vdv*e_i+u*2*3*z"2*d_iz*vdv*e_i+u*2*z"3*d_ivdv*e_i+u*2*z*3*v~ 
dv*d_ie_i 

uu2:=-2*u*d_iu*z*3*e_i-u*2*3*z*2*d_iz*e_i-u*2*z"3*d_ie_i 

uu3 : =+2*u*d_iu*z~2*v_i+u~2*2*z*d_iz*v_i+u~2*z~2*d_iv_i 

uu4:=-d_iu*z~3*eda*e_i-u*3*z~2*d_iz*eda*e_i-u*z"3*d_ieda*e_i-u*z~3*eda*d_ie_i 
II II 

"Thus second peirtial (D_i)~2 of u*z is" 
dd_iphi : =uul+uu2+uu3+uu4 

"Expanding right side yields" 

dd_iphi:=3*u"3*z'5*vdv'2*e_i"2-6*u"3*z"5*vdv*e_i"2+3*u"3*z"5*e_i'2-u"3*z"4*ee_~ 
ii*vdv+u*3*z*4*ee_ii+u*3*z*4*ev_ii*vdv-u"3*z*4*ev_ii+u*3*z"4*vdv*e_i"2+5*u*3*z~ 
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~4*vdv*e_i*v_i-u~3*z~4*e_i~2-5*u~3*z~4*e_i*v_i+u~3*z~3*un_ii*vdv-u~3*z~3*un_ii~ 

+2*u~3*z~3*v_i~2-6*u~2*z~5*eda*vdv*e_i"2+6*u"2*z"5*eda*e_i"2+u'2*z"4*eda*ee_ii~ 

-u~2*z~4*eda*ev_ii-u~2*z~4*eda*e_i*2-5*u*2*z~4*eda*e_i*v_i-3*u~2*z"4*vda*e_i"2~ 

-u*2*z*3*eda*im_ii-u*2*z*3*ea_ii-u"2*z*3*e_i*a_i+3*u*z*5*eda"2*e_i"2+u*z*4*edd~ 

ot_a*e_i*2 

II II 

"After summation over index i we get the Laplacian of Phi=u*z" 

lp_phi:=+3*u"3*z'5*vdv"2*ede-6*u"3*z~5*vdv*ede+3*u"3*z'5*ede-u"3*z"4*ede*vdv+u~ 

-3*2;-4*ede+u*3*z*4*edv*vdv-u*3*z"4*edv+u"3*z*4*vdv*ede+5*u"3*z*4*vdv*edv-u*3*z~ 

*4*ede-5*u*3*z*4*edv+u*3*z*3*3*vdv-u*3*z*3*3+2*u*3*z*3*vdv-6*u*2*z*5*eda*vdv*e~ 

de+6*u^2*z"5*eda*ede+u'2*z'4*eda*ede-Ti"2*z"4*eda*edv-u'2*z"4*eda*ede-5*u'2*z"4~ 

*eda*edv-3*u'2*z"4*vda*ede-u*2*z"3*eda*3-u"2*z*3*eda-u*2*z"3*eda+3*u*z*5*eda"2~ 

*ede+u*z"4*eddot_a*ede 
II II 

"Check if wave equation is satisfied" 
wave : =lp_phi-dd_phi 
"Expanding the right side yields" 
wave : =0 

"Eureka! It worked." 
II II 

"For the record expaoid the expression defining variable lp_phi" 
lp_phi:=3*u"3*z~5*vdv"2-6*u"3*z*5*vdv+3*u*3*z"5+6*u*3*z"4*vdv-6*u"3*z*4-u"3*z*~ 
3*vdv+3*u"3*z"3-6*u"2*z"5*eda*vdv+6*u"2*z*5*eda-6*u*2*z"4*eda-3*u"2*z*4*vda+u*~ 
2*z"3*eda+3*u*z"5*eda"2+u*z"4*eddot_a 

Check that vector potential A = uzv satisfies the wave equation 

"Check that vector potential aa=u*z*v satisfies wave equation" 
11 11 

"Generate Laplacian of vector potential aa=u*z*v " 
"Compute first partials D_i_aa=D_i [uzv] " 
"d_i_aa= d_iu*z+v +u*d_iz*v +u*z*d_iv" 
"In expanded form the expression has 5 terms" 

"d_i_aa= -z~3*e_i*eda*u*v +z~3*e_i*u~2*v*vdv -z*3*e_i*u"2*v ~ " 
" -z"2*a*e_i*u +z"2*u"2*v*v_i" 

"Compute D_i derivative of each term of the above sum" 

ttl :=-3*z~2*d_iz*e_i*eda*u*v-z~3*d_ie_i*eda*u*v-z~3*e_i*d_ieda*u*v-z~3*e_i*eda~ 
*d_iu*v-z~3*e_i*eda*u*d_iv 

tt2:=+3*z"2*d_iz*e_i*u"2*v*vdv+z*3*d_ie_i*u"2*v*vdv+z"3*e_i*2*u*d_iu*v*vdv+z"3~ 
*e_i*u"2*d_iv*vdv+z"3*e_i*u"2*v*d_ivdv 

tt3:=-3*z'2*d_iz*e_i*u"2*v-z~3*d_ie_i*u"2*v-z"3*e_i*2*u*d_iu*v-z"3*e_i*u*2*d_iv 

tt4:=-2*z*d_iz*a*e_i*u-z"2*d_ia*e_i*u-z"2*a*d_ie_i*u-z~2*a*e_i*d_iu 

tt5 : =+2*z*d_iz*u~2*v*v_i+z~2*2*u*d_iu*v*v_i+z"2*u*2*d_iv*v_i+z*2*u"2*v*d_iv_i 
II II 

d_i_d_i_aa:=ttl+tt2+tt3+tt4+tt5 
II II 

"Expanded the right side of above expression yields" 

d_i_d_i_aa: =3*z*5*e_i*2*eda*2*u*v-6*z*5*e_i*2*eda*u*2*v*vdv+6*z"5*e_i*2*eda*u*" 
2*v+3*z*5*e_i*2*u*3*v*vdv*2-6*z*5*e_i*2*u*3*v*vdv+3*z"5*e_i*2*u"3*v+3*z"4*a*e_~ 
i~2*eda*u-3*z~4*a*e_i"2*u~2*vdv+3*z~4*a*e_i~2*u"2-z*4*e_i~2*eda*u"2*v+z"4*e_i*~ 

2*eddot_a*u*v+z~4*e_i~2*u~3*v*vdv-z~4*e_i~2*u~3*v-3*z~4*e_i~2*u~2*v*vda-5*z"4*~ 
e_i*eda*u*2*v*v_i+5*z*4*e_i*u*3*v*v_i*vdv-5*z*4*e_i*u*3*v*v_i+z*4*eda*ee_ii*u*"' 
2*v-z"4*eda*ev_ii*u"2*v-z"4*ee_ii*u"3*v*vdv+z"4*ee_ii*u"3*v+z*4*ev_ii*u"3*v*vd~ 
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v-z~4*ev_ii*u~3*v-z~3*a*e_i~2*u~2-3*z~3*a*e_i*u~2*v_i+z~3*a*ee_ii*u'2-z~3*a*ev~ 
_ii*u~2-z'3*a_i*e_i*u'2*v+z~3*dot_a*e_i~2*u-z"3*ea_ii*u~2*v-z"3*eda*u"2*im_ii*~ 

v+z*3*u*3*im_ii*v*vclv-z*3*u*3*iiii_ii*v+2*z"3*u"3*v*v_i~2-z~2*a*u~2*un_ii 
ti It 

"To compute Laplaciaii(aa) we need to take sum over index i" 

"remembering that sum(e_i*e_i)=ede; sum(e_i*v_i)=edv; sum(ee_ii)=ede ; " 

"sum(un_ii)=3; s\m(v_i*2)=vdv etc." 
II II 

"To this end edit previous expression to make proper changes" 

"Thus Laplacian(aa) is" 

lp_aa : =3*z"5*ede*eda"2*u*v-6*z"5*ede*eda*u*2*v*vdv+6*z"5*ede*eda*u*2*v+3*z*5*'' 

ede*u"3*v*vdv'2-6*z~5*ede*u"3*v*vdv+3*z~5*ede*u~3*v+3*z~4*a*ede*eda*u-3*z~4*a*~ 

ede*u*2*vdv+3*z*4*a*ede*u*2-z*4*ede*eda*u"2*v+z"4*ede*eddot_a*u*v+z~4*ede*u"3*~ 

v*vdv-z~4*ede+u~3+v-3*z~4*ede*u~2*v*vda-5*z*4*edv*eda+u"2*v+5*z"4*edv*u"3*v*vd~ 

v-5*z~4*edv*u~3*v+z~4*eda*ede*u~2*v-z~4*eda*edv*u~2*v-z~4*ede*u~3*v*vdv+z~4*ed~ 

e*u*3*v+z*4*edv*u*3*v*vdv-z*4*edv*u*3*v-z*3*a*ede*u~2-3*z~3*a*edv*u*2+z*3*a*ed"' 

e*u*2-z*3*a*edv*u~2-z~3*eda*u*2*v+z"3*dot_a*ede*u-z~3*eda*u*2*v-z"3*eda*u"2*3*~ 

v+z*3*u"3*3*v*vdv-z"3*u"3*3*v+2*z"3*u"3*v*vdv-z"2*a*u"2*3 

n II 

"Now compute time derivative D_aa of [u*z*v] " 
d_aa : =+d_u*z*v+u*d_z*v+u*z*d_v 
"It expands to" 

d_aa: =+u~2*z~3*v* (1-vdv) -u"2*z"2*v+eda*u*z"3*v+a*u*z"2 
11 II 

"Compute derivative of each term of above expression" 

ppl:=+2*u*d_u*z"3*v*(l-vdv)+u*2*3*z"2*d_z*v*(l-vdv)+u"2*z'3*d_v*(l-vdv)-u"2*z"~ 
3*v*d_vdv 

pp2 :=-2*u*d_u*z"2*v-u~2*2*z*d_z*v-u~2*z~2*d_v 

pp3 : =+d_eda*u+z*3+v+eda*d_u*z*3*v+eda*u*3*z*2*d_z*v+eda*u*z"3*d_v 

pp4 : =+d_a*u*z"2+a*d_u*z"2+a*u*2*z*d_z 
II II 

"The second derivative of aa with respect to time is" 

dd_aa : =ppl+pp2+pp3+pp4 
II II 

"It expands to" 

dd_aa:=3*u"3*z"5*vdv"2*v-6*u"3*z*5*vdv*v+3*u*3*z"5*v+6*u~3*z"4*vdv*v-6*u*3*z"4~ 
*v-u"3*z"3*vdv*v+3*u"3*z"3*v-6*u"2*z*5*vdv+eda*v+6*u"2*z*5*eda*v-3*u"2*z*4*vdv~ 

*a-6*u'2*z~4*eda*v-3*u~2*z~4*vda*v+3*u"2*z~4*a+u~2*z~3*eda*v-4*u*2*z"3*a+u*2*z~ 

~2*a+3*u*z"5*eda"2*v+3*u*z"4*eda*a+u*z~4*v*eddot_a+u*z'3*dot_a 
11 II 

"D' Alembertian of aa is" 
wave_aa : =lp_aa-dd_aa 

"Expanding right side of above yields:" 

wave_aa:=0 
11 II 

"Thus the wave equation is satisfied!" 
11 II 

"For the record expand the right side of statement defining variable dd_aa" 

dd_aa:=3*u*3*z"5*vdv*2*v-6*u"3*z~5*vdv*v+3*u"3*z"5*v+6*u"3*z"4*vdv*v-6*u"3*z*4~ 

*v-u*3*z~3*vdv*v+3*u~3*z~3*v-6*u*2*z*5*eda*vdv*v+6*u"2*z*5*eda*v-6*u"2*z*4*eda~ 

*v-3*u"2*z~4*vda*v-3*u~2*z~4*vdv*a+3*u'2*z~4*a+u~2*z"3*eda*v-4*u~2*z"3*a+u~2*z~ 

~2*a+3*u*z~5*eda~2*v+3*u*z~4*eda*a+u*z"4*eddot_a*v+u*z"3*dot_a 
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Check that Lorentz gauge formula V • A + D(j) = is satisfied. 

"Check Lorentz gauge formula Div(A)+D(Phi)=0" 
II II 

"Computing D_iA_i=D_i [uzv_i] " 

d_i_a_i : =d_iu*z*v_i+u*d_iz*v_i+u*z*d_iv_i 
II II 

"Expanding the right side of above statement yields" 

d_i_a_i :=u*2*2*3*vdv*e_i*v_i-u*2*z"3*e_i*v_i+u*2*z"2*v_i*2-u*z*3*eda*e_i*v_i-u~ 

*z*2*ea_ii 
II II 

"Thus the divergence of the vector field A is" 

div_a:=u*2*z*3*vdv*edv-u*2*z*3*edv+u*2*z"2*vdv-u*z"3*eda*edv-u*z*2*eda 
II II 

"Time derivative of Phi is" 

d_phi : =d_u*z+u*d_z 
II II 

"Lorentz gauge is" 

gauge : =div_a+d_phi 
II II 

"Expanding the right side of above yields" 

gauge :=0 
II II 

"For the record expand the expression defining variable d_phi" 
d_phi:=-u"2*z"3*vdv+u"2*z"3-u"2*z*2+u*z"3*eda 

Compare Feynman's field with Lienard-Wiechert's field 

"Compare Feynman's electric field f f =u"2*e+u" (-l)*D(u'2*e)+D"2e" 

"with Lieneird-Wiechert electric field ee=-grad(u*z)-D(u*z*v) " 
II II 

"First term of the sum defining ff is" 

ttl :=u~2*e 

"The second term is" 

tt2 : =u* (-1) * (2*u*d_u*e+u*2*d_e) 

"To compute the third term notice first that" 

d_e : =u*z* (e-v) -e*u 

"Thus for the third term we have the formula" 

tt3 : =d_u*z* (e-v) +u*d_z* (e-v) +u*z* (d_e-d_v) -d_e*u-e*d_u 
II II 

"Thus we have for Feynman field ff" 

ff :=ttl+tt2+tt3 
II II 

"Now compute components ee_i of the electric field ee" 
"where ee=-grad(u*z)-D(u*z*v) " 

ee_i : =- (d_iu*z+u*d_iz) - (d_u*z*v_i+u*d_z*v_i+u*z*d_v_i) 
"The above expamds to" 

ee_i :=-u~2*z*3*vdv*e_i+u*2*z*3*vdv*v_i+u"2*z*3*e_i-u"2*z*3*v_i+u*z"3*eda*e_i-u~ 

*z*3*eda*v_i-u*z*2*a_i 
II II 

"Thus in vector notation ee is" 

ee :=-u'2*z'3*vdv*e+u~2*z"3*vdv*v+u~2*z~3*e-u"2*z~3*v+u*z"3*eda*e-u*z~3*eda*v-u~ 
*z*2*a 
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difference : =f f -ee 

"Expanding in the above statement the expression on the right side yields" 
difference : =0 

"So the electric fields coincide!" 
It II 

"Formulas for magnetic fields involve cross product, programming of which is" 
"cumbersome. It is easier to do it by a direct computation as in the text." 
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